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QUANTUM ERGODICITY OF BOUNDARY VALUES OF 
EIGENFUNCTIONS 

ANDREW HASSELL AND STEVE ZELDITCH 


Abstract. Suppose that Q C is a bounded, piecewise smooth domain. 
We prove that the boundary values (Cauchy data) of eigenfunctions of the 
Laplacian on H with various boundary conditions are quantum ergodic if the 
classical billiard map /3 on the ball bundle B*(dQ) is ergodic. 

Our proof is based on the classical observation that the boundary values of 
an interior eigenfunction (pxi is an eigenfunction of an operator 

Fh on the boundary of Q with h = . In the case of the Neumann bound¬ 

ary condition, is the boundary integral operator induced by the double 
layer potential. We show that F^ is a semiclassical Fourier integral operator 
quantizing the billiard map plus a ‘small’ remainder; the quantum dyanmics 
defined by F^ can be exploited on the boundary much as the quantum dynam¬ 
ics generated by the wave group were exploited in the interior of domains with 
corners and ergodic billiards in the work of Zelditch-Zworski (1996). Novel¬ 
ties include the facts that F^ is not unitary and (consequently) the boundary 
values are equidistributed by measures which are not invariant under (3 and 
which depend on the boundary conditions. 

Ergodicity of boundary values of eigenfunctions on domains with ergodic 
billiards was conjectured by S. Ozawa (1993), and was almost simultaneously 
proved by Gerard-Leichtnam (1993) in the case of convex domains (with 
continuous tangent planes) and with Dirichlet boundary conditions. Our meth¬ 
ods seem to be quite different. Motivation to study piecewise smooth domains 
comes from the fact that almost all known ergodic domains are of this form. 
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1. Introduction 

The purpose of this article is to prove quantum ergodicity 

(1.1) (AhjUjjUj) / o'(A) dfiB, j ^ oo along a density one sequence 

J B*Y 

of the boundary values of interior eigenfunctions 

{ A Uj = A| Uj in n, {uj,Uk)L^n) = Sjk, 

Buj\Y = 0, Y = dn 

of the Euclidean Laplacian A on a compact piecewise smooth domain C K." and 
with classically ergodic billiard map f3 : B*Y B*Y, where Y = dfl. Here Ah is 
a zeroth order semiclassical pseudodifferential operator on Y. The relevant notion 
of boundary values (i.e. Cauchy data) Uj depends on the boundary condition B, 
as does the classical limit measure dfiB according to which the boundary values 
become equidistributed. Our methods cover Dirichlet, Neumann, Robin and more 
general boundary conditions of the form 

(1.2) Bu = d^u — K{u\y), if S il'^(F) self-adjoint with non-negative symbol 

which we refer to as a d'^-Robin boundary condition. 

For the interior problem, ergodicity of eigenfunctions of Laplacians on bounded 
domains with corners and with ergodic billiard flow was proved by Zelditch-Zworski 
|ZZw| . following an earlier proof by Gerard-Leichtnam El] in the case of 
convex domains. Our proof of boundary ergodicity is independent of these proofs 
in the interior case. In the case of manifolds without boundary, results on ergodicity 
of eigenfunctions originate in the work of A. I. Schnirelman |Sch| and were carried 
forward by Zelditch | 2 ] and Colin de Verdiere ITWI . We refer to |Z2| for background 
and a simple proof which will be developed here. 

Ergodicity of boundary values of eigenfunctions was conjectured by S. Ozawa 
im in 1993 and was independently stated and proved by Gerard-Leichtnam ELI in 
the same year in the case of Dirichlet boundary conditions on convex domains 
(i.e. domains whose unit normal is Lipschitz regular). Our extension to nonconvex 
piecewise smooth domains is motivated by the fact that ‘most’ known ergodic bil¬ 
liard tables are nonsmooth. More precisely, convex domains never have ergodic 
billiards, and the only known convex ergodic billiard tables appear to be the 
Bunimovich stadium example fBU^[BU3| and its relatives in two dimensions (i.e. 
plane domains formed by straight segments and circular segments). In particu¬ 
lar, no convex ergodic domains are known in dimensions > 3 (the example 
cited in m from |BU3| is not ergodic) and there is some doubt that they exist 
Em]. The known higher dimensional ergodic stadia (i.e. with only focussing or 
nowhere dispersing boundary faces, consisting of convex or flat components) are 
non-convex and have corners [BR.1 IBR2| . A sizable collection of ergodic, piecewise 
smooth plane domains is furnished by generic polygons EMS]. Many other exam¬ 
ples of non-convex ergodic (and hyperbolic) billiard domains in higher dimensions 
are given by dispersing billiard domains bounded by unions of concave boundary 
components (see e.g. m)- 

Our proof is based on a reduction to the boundary of the eigenvalue problem. 
The intuitive idea is that the Cauchy data (wjlv, di,Uj\Y) of interior eigenfunctions 
Uj provide a kind of quantum cross section to the interior eigenfunctions, just as 
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the billiard map P on B*Y provides a kind of cross section to the billiard flow on 
T*n. To be precise, our starting point is the classical observation that the boundary 
value of an interior eigenfunction with eigenvalue is an eigenfunction of a certain 
boundary operator Fh with h = A“^. We analyse Fh (on a convex domain, or a 
modification of Fh for a nonconvex domain) as a semiclassical Fourier integral oper¬ 
ator quantizing /3 plus a remainder which is almost ignorable. Boundary ergodicity 
has some new features which are not present in interior ergodicity, stemming from 
the fact that the dynamics generated by Fh defines an endomorphism but not an 
automorphism of the observable algebra. 

To state our results, we will need some notation: Let Ag denote the positive 
Laplacian on n with boundary conditions Bu = 0. Then has discrete spec¬ 
trum 0 <Ai<A 2 <---^oo, where we repeat each eigenvalue according to its 
multiplicity, and for each \j we may choose an normalized eigenfunction Uj . 

To each boundary condition B corresponds 

• A specific notion of boundary value of the eigenfunctions Uj. We denote 
the L^-normalized boundary values by u'j = u^j/\\u^\\. 

• A specific measure d/rs on B*{Y). 

• A specific state ojb on the space di)))!") of semiclassical pseuodifferential 
operators of order zero. 

The correspondence is dictated by the local Weyl law ('T;emma |l .211 for the bound¬ 
ary condition B. Here is a table of the relevant boundary value notions. In the 
table, K, denotes a (7°° function on Y while k is the principal symbol of the operator 
K G Il'^(y) in 11.211 . and da is the natural symplectic volume measure on B*Y. We 
also define the function 7 (g) on B*Y by 

(1-3) y{q) = Vl - hp, q={y,r]). 


Boundary Values 

B 

Bu 


dyLB 

Dirichlet 

u\y 

dyU\Y 

j{q)da 

Neumann 

duU\Y 

u\y 

y{q)-'-da 

Robin 

{dyU — ku)\y 

u\y 

b 

1 

'F^-Robin 

{dyu - Ku)\y 

u\y 

y{q)da 

yiqY + k{qY 


Throughout, we assume our domain H is a piecewise smooth manifold embedded 
in R". Hence its boundary, denoted Y, is the union of a smooth part Y° and a 
singular set E, which has measure zero (see section|3for background on such mani¬ 
folds). The metric on H is understood to be the Euclidean metric. Our main result 
is that, if the billiard ball map (3 on B*Y is ergodic, then the boundary values of 
eigenfunctions are quantum ergodic. As reviewed in § 2 , quantum ergodicity has to 
do with time and space averages of observables. The relevant algebra of observables 
in our setting is the algebra of zeroth order semiclassical pseudodifferential 

operators on Y, depending on the parameter h G [0, ho]. We denote the symbol of 
A = Ah G 'I')](Y) by a = a{y,r],h). Thus a{y,ri) = a{y,ri,0) is a smooth function 
on T*Y. We further define states on the algebra 'I'°(Y) by 


ujb{A) = 


vol(5'" ^)vol(H) 


a{y,v)df^B- 


(1.4) 


B*Y 
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Then our main result is 


Theorem 1.1. Let LI C K" be a bounded piecewise smooth manifold (see Def¬ 
inition ro) with ergodic billiard map. Let {mJ} be the boundary values of the 
eigenfunctions {uj} of Ab on L‘^{Ll) in the sense of the table above. Let Ah be 
a semiclassical operator of order zero on Y. Then there is a subset S of the positive 
integers, of density one, such that 


(1.5) 


lim {Ah.jU^.j,u^f) = ujb{A), B= Neumann, Robin or -Robin, 
j^oo,jeS ■' 

lim \~'^{Ah.u\,u^A = ujb{A), B= Dirichlet, 

j^oojGS J J 


where hj = Xj ^ and lob is as in (US- 


• Let us give the results more explicitly for the identity operator. For the Neumann 
boundary condition, 


lim I 

j—too.jGS 


Il2(y) 


2 voi(y) 
voi(n) ’ 


while for the Dirichlet boundary condition, 


lim A, 




2voi(r) 

n vol(D) 


• As mentioned above, the case of Dirichlet boundary conditions on convex 
domains was proved earlier by Gerard-Leichtnam m Their proof was based on 
an identity (En, Theorem 2.3) relating quantum limit measures of interior Dirich¬ 
let eigenfunctions to those of its boundary values. As they point out, the proof 
assumes regularity and does not apply to domains with corners. Subsequent 
to the initial version of this article (which only proved ergodicity for convex billiard 
domains), N. Burq m proved boundary ergodicity for all Riemannian manifolds 
with corners (including non-convex domains) and all boundary conditions consid¬ 
ered here. His method is to extend the method of to general piecewise smooth 
domains and boundary conditions, and thus to reduce the proof of boundary quan¬ 
tum ergodicity to the known interior quantum ergodicity result in |ZZw| . His proof 
also uses the results of the present article on non-concentration of eigenfunctions 
at the corners. 

• The Neumann and Dirichlet limit measures dfj.B can be understood as follows: 
First, 7 “^dcr is the projection to B*Lt of the Liouville measure on the set 5'*„y 
of inward pointing unit vectors to Ll along Y under the projection taking a vector 
to its tangential component. We also recall that (3 is symplectic with respect to 
the canonical symplectic form da on B*Y. Thus, boundary values of Neumann 
eigenfunctions are equidistributed according to the measure on B*Y induced by 
interior Liouville measure rather than the boundary symplectic volume measure. 
In the case of Dirichlet boundary conditions, the boundary value u’j is taken to 
be normal derivative of the eigenfunction at the boundary. The symbol of hd,y, 
restricted to the spherical normal bundle, and then projected to B*Y is equal to 
7 , so we should expect to get the square of this factor in the Dirichlet case (since 
mu is quadratic in uj) compared to the Neumann case. The normal derivatives 
in the Dirichlet case also account for the factor in 
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The fact that the quantum limit state ojb and the corresponding measure 
do not in general coincide with the natural symplectic volume measure da on B*Y 
will be traced in the fact that the quantum dynamics is defined by an endo¬ 
morphism rather than an automorphism of the observable algebra. Let us explain 
how this works in the case of Neumann boundary conditions. In this case, the 
dynamics are generated by the operator Fh on Y with kernel 

Fhiy.y') = ‘2-^Go{y,y',h~^), y:^y'eY, where 

( 1 . 6 ) , 

Go{y, y', A) = - z'|)-(-2)/2 HaW^_,(A|z - z'|) 

is the free outgoing Green function on R”. By virtue of Green’s formula 

(1.7) Uj{z) = J [d^^,Go{z, y', Xj)uj{y') - Go{z, y', Xj)d^^,Uj{y')] da{y') 
for any solution of IXuj = X'juj, and the jump formula 

( 1 . 8 ) limv- / ‘^duy'Go{z,y',Xj)uj{y')da{y') = Uj{y) + Fh{uj){y), 

z^y&Y Jy “ 

this operator leaves the boundary values of Neumann eigenfunctions invariant: 

(1.9) Fh,u) = u), j = 1,2,... 

It follows that the states 

(1.10) p,{A) ■■= 

are invariant for Fh^ ■ Similar invariance properties hold for the other boundary 
conditions. As we will show, the family {F/^} defines a semiclassical Fourier integral 
operator associated to the billiard map /3 (for convex 17), plus some terms which 
turn out to be negligible for our problem. The quantum dynamics on is thus 

generated by the conjugation 

( 1 . 11 ) = Fh, Ah, Fn^ 

This is analogous to the interior dynamics generated by 

(1.12) at{A) = UtAU:, U{t) = , 

but it has one important difference: unlike U(t ), Fh is not unitary or even normal. 
Indeed, the zeroth order part of F^Fh is a semiclassical pseudodifferential operator 
with a non-constant symbol. This new feature of the quantum ergodicity problem 
is one of the prinicipal themes of the present article. 

We now outline the proof, emphasizing the aspects which are new to the bound¬ 
ary case. The general strategy is the same as in and relies on two main 

ingredients: a local Weyl law for the u^, and an Egorov type theorem for an operator 
Fh- Naturally, the reduction to the boundary brings in additional considerations, 
which are of some independent interest. 

We begin with the local Weyl law, which has nothing to do with ergodicity; it is 
valid for all domains Q. 
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Lemma 1.2. Let be either the identity operator on Y or a zeroth order semi- 
elassieal operator on Y with kernel supported away from the singular set. Then for 
any of the above boundary conditions B, we have: 

(1.13) 


lim 

X^oo 


1 


A,<A 


lob(A), B= Neumann, Robin or -Robin, 


lim 
A—^oo 


1 


A,<A 


B = Dirichlet. 


When Al is a multiplication operator, and for Dirichlet boundary conditions, this 
local Weyl law was essentially proved by Ozawa El. In §3 we extend the proof to 
general semiclassical pseudodifferential operators and to the boundary conditions 
described above. That allows us to capture uniform distribution of eigenfunctions 
in phase space rather than just in configuration space. For multiplication operators, 
the local Weyl law can be obtained from the interior local Weyl law by Hadamard’s 
variational formula with respect to the boundary conditions. This observation was 
first made by Ozawa f0^[(^l04| . 

From the local Weyl law we deduce an invariance property of the limit states. For 
notational simplicity, we confine ourselves here to Neumann boundary conditions, 
where the boundary operator is CS; analogous invariance properties hold for other 
boundary conditions with small modifications to F^. 


Corollary 1.3. The state WNeu is invariant under Fh-' ujf^eu{F^AFh) = a;Neu(^). 

Indeed, the states Pj{A) = {Ah^UpUj) are invariant so any average or limit of 
averages of these states will be invariant. 

The Egorov type result for the operator Ff, is as follows: 


Lemma 1.4. Let D C K” be a bounded piecewise smooth convex domain, let /3 
denote the billiard map on B*Y° and let Ah = Op(a/i) be a zeroth order operator 
whose symbol a{y, 77, 0 ) at h = 0 is supported away from 

{|r?| = l}UEur'(S). 

Let 7 be given by (11.311 . Then 

F^AhFh = Ah+ Sh, 


where Ah is a zeroth order pseudodifferential operator and ||S';i||L 2 ^i 2 < Ch. The 
symbol of Ah is 


(1.14) 


0 , 


q G B*Y 
q i B*Y. 


This is a rigorous version of the statement that Fh quantizes the billiard ball 
map. The unusual transformation law of the symbol reflects the fact that cm is 
not an automorphism. This Egorov theorem is relevant to the Neumann boundary 
problem. In the Dirichlet case, the relevant operator is In the Robin case there 
is a lower order term, while in the 'h^-Robin case there is a second term of the same 
order. For nonconvex domains, we replace Fh by a modified invariant operator 
which removes the spurious wavefront set of Fh not associated with fd. 

We now sketch the completion of the proof in the case of the Neumann boundary 
condition. As in the case of automorphisms, it is essentially a convexity argument 
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(see §2 or |Z2| for this point of view). For simplicity of exposition, we temporarily 
ignore the problems caused by the corners and pretend that the domain is smooth; 
details on the corner issues appear in Section [7| 

To show that 


{Aulu^^ 


^jl ' ^Neu(-^); 

along a density one subsequence of integers j is essentially to show that 

1^1 . . 2 


(1.15) 


limsup 
A^oo N{\) 




((A - WNeu(^))w' 




= 0 . 


Due to the novel form of the local Weyl law and the Egorov theorem, we cannot 
apply the averaging argument of |Z2| directly to H1.15II . Instead, we introduce an 
auxiliary pseudodifferential operator R which is almost invariant under conjugation 
by Fh- It suffices that its symbol is given (approximately) by 

vol(S'”-i)vol(D) 


(1.16) 


o-R(g) C 7 (q) = c(l - |77l^)^^^ 


Avo\{B*Y) 

As an intermediate step, we prove the following analogue of 11.1511 : 


Lemma 1.5. Let LI C ffi." be a bounded piecewise smooth domain, with ergodic 
billiard map, and let Ah be a zeroth order operator. For all e > 0, there exists a 
pseudodifferential operator Rh of the form U.lbV such that 


(1.17) 


lim sup 

A—^oo 


1 

W) 


E 

A,<A 


WNeu(^) Rhj)Upu’’j) 


2 

< e. 


The proof uses Cauchy-Schwarz to bound the left side of da by 

1 


(1.18) 


limsup 
A^oo 7V(A) 


^(((A-i?C.Neu(A))^w5,u5)). 


A,<A 


We then use the invariance properties and Corollary to replace A by the 
time-average A^ defined by 


1 ^ 

{AN)H = j^Y.^{FtrAhFt 

^ k=l 

After replacing A by A^v, it follows from the local Weyl law of Lemma ^21 that the 
limit in is given by 

(1.19) WNeu((AAr — WNeu(A)i?)^) = y ^Cr(AAr) — lA;Neu(A)CT(i?)^ dliNeu- 

By applying Lemma lL^ iteratively we see that the the symbol cr(AAr) of Ajv is 7 (g) 
times a discrete average of 7 “^ct(A) over iterates of the billiard map. Hence, this 
converges to 7 (g) x f 7 “^cr(A) dcr/ f 1 da, which is equal to C 7 ((;)a;Neu(A). This is 
approximately equal to the symbol of a;Neu(A)i?. Thus, the integral 11.1911 becomes 
small as N 00 , and we can make arbitarily small by choosing a sufficiently 

good operator R. We deduce 11.151) by considering (tTTTll with A replaced by 
WNeu(A)Id, and subtracting this from 11.1711 . 

In this sketch, we implicitly assumed that Egorov’s theorem and other pseudodif¬ 
ferential manipulations were valid on all of Y. However, they are only valid away 
from the corners and hence we have only obtained an equidistribution law away 
from the corners. But Lemmawith Ah = Id, shows that eigenfunctions do not 
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concentrate at the corners and therefore WNeu is the full limit measure. The proof 
is similar in the case of the other boundary conditions. 

Let us now describe the organization of this article. Although the details of the 
arguments differ somewhat (and sometimes significantly) as the boundary condition 
B varies, we present complete details first in the case of the Neumann boundary 
condition, and then describe the necessary modifications for the other boundary 
conditions. Robin boundary condition comes second because it is a deformation of 
the Neumann boundary condition and the details are quite similar. We then take 
up the Dirichlet boundary condition, which is singular relative to Neumann. The 
dt^-Robin boundary condition provide a bridge between Neumann and Dirichlet and 
require the most serious modifications; we present the details for these conditions 
last. We shall also assume that our domain is convex until the Section CH where 
we adapt the argument to nonconvex domains. To guide the proof, we also start 
in Section |2 with an overview of how boundary ergodicity differs in principle from 
the general structure of proof of quantum ergodicity given in m In the appendix 
we summarize properties of heat kernels needed for the proof. 

We wish to thank Ben Andrews, Nicolas Burq, Rafe Mazzeo, and Tom ter Elst 
for helpful discussions; Leonid Bunimovich and Maciej Wojtkowski for informing 
us about the state of the art on Euclidean domains with ergodic billiards; Maciej 
Zworski for advising us that the generalization from convex to general domains 
should not be hard and encouraging us to treat this case; and Alan McIntosh, 
Andreas Axelsson, Monique Dauge and Michael Taylor for background on Lipschitz 
domains. 


2. Quantum ergodicity of endomorphisms 

Quantum ergodicity is concerned with quantizations of classically ergodic Hamil¬ 
tonian systems. The typical example is the wave group Ut = of a compact Rie- 
mannian manifold (M, g) without boundary, whose geodesic flow G* : S*M ^ S*M 
is ergodic on the unit co-sphere bundle with respect to Liouville measure dg,L- Er¬ 
godicity on the quantum level is a relation between the time average 

(^):=w- lim ^ ( UtAUtdt 
1 —^oo Al. J _ i-p 

of an observable A S iI'°(M) and the space average^ defined to be the constant 
operator 

uj{A) I, where w(A) = . / crAd^iL■ 

VOl{S*M) Jg,M 

Here, iI'°(M) is the space of zeroth order pseudodifferential operators over M and 
a A is the principal symbol of A. 

The system is said to be quantum ergodic if 

(2.1) (A)= u;{A)I + K, with ^ 0, 

where || • \ \hs is the Hilbert-Schmidt norm and where H^ is the spectral projection 
onto the span of eigenfunctions of -s/A of eigenvalue < A. Thus, time and space 
averages agree up to an operator K which is negligible in the semi-classical limit. 
By expressing all operators in terms of the eigenfunctions of t/j, one sees that this 
formulation is equivalent to 11.1511 . 




QUANTUM ERGODICITY OF BOUNDARY VALUES OF EIGENFUNCTIONS 


9 


The quantum dynamics is thus the C* dynamical system defined by 

the one-parameter group of automorphisms 

( 2 . 2 ) at{A)=UtAUl 

on the norm closure of iho- It is a general fact that (K, at, iIi®) is quantum ergodic as 
long as u! is an ergodic state |Z2| . After applying the Schwarz inequality for states 
and the Egorov theorem, the proof is reduced to the mean ergodic theorem on 
the symbol level. 

This general fact does not apply immediately to the wave group Ut of the Lapla- 
cian As on a bounded domain C M in a Riemannian manifold {M,g), with one 
of the above boundary conditions B. The problem is that conjugation by the wave 
group does not quite define an automorphism of any natural algebra of pseudodif¬ 
ferential operators. Nevertheless, by using suitable approximations, one can prove 
that n~T^ holds for pseudodifferential operators A on the ambient manifold which 
are essentially supported in the interior of |ZZw| . 

The situation we consider in this paper is complicated in one further way: As 
mentioned in the introduction, the relevant quantum dynamics cm defines an 
endomorphism rather than an automorphism of the relevant algebra of observables. 
Let us first explain this on the classical level. For simplicity, we restrict to smooth 
plane domains and Neumann boundary conditions. The general case is similar. 

Let n C K" be a smooth domain and let 


(2.3) [3 : B*Y B*Y 

denote the billiard map (see Sectional . In the spectral theory of dynamical systems, 
one studies the dynamics of [3 through the associated ‘Koopman’ operator 

U : L\B*Y, da) ^ L\B*Y, da), Uf{0 = /(/3(C)). 


Here, da denotes the usual /3-invariant symplectic volume measure on B*Y. From 
the invariance it follows that lA is a. unitary operator. When /3 is ergodic, the unique 
invariant L^-normalized eigenfunction is a constant c, and one has the mean ergodic 
theorem (see e.g. |P], chapter 2 ) 


(2.4) 


1 


N 

lim Wi-^W 
N—*oo N 

n=l 


(/)-(/,c)||- 0 . 


Suppose now that we introduce a positive function 7 G C°°{B*Y), and define a 
new operator T by 

Then T is not unitary on L‘^{B*Y,da). This is the situation we find ourselves in 
by virtue of the Egorov Lemma d in the case of Neumann boundary conditions. 
We find that the orthogonal projection P onto the invariant functions for T 
relative to the invariant inner product has a somewhat different form: 


Proposition 2.1. Consider the operator T on L^{B*Y,da). We claim: 

• (i) The unique positive T-invariant function is given by 7 . 

• (a) T is unitary relative to the inner product {{,)) on B*Y defined by the 
measure dv = ^~‘^da. 
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• (Hi) When (3 is ergodic, the orthogonal projection P onto the invariant 
L'^-eigenvectors has the form 


P{f) 


((/>7)),,_ 1 

(( 7 , 7 ))^ vol(B*y) 


where c is as in OHIl . 



/7 ^da] 7 = cwNeu(/) 7 


Proof. An measure pda is invariant under T if and only if 

7ifL d(/3(C)) = p{C) 7"^(C)p(C) is invariant under (3. 

The unique positive solution (and the unique solution in the ergodic case) is given 
by pda = 'yda. This proves (i). 

To prove (ii), we seek an invariant inner product ((,)) of the form 


((/. 5 ))= [ fiOgiOdv- 

JB 

We note that 

= 0 / 3-1 

= {{f,T-^g)) ^ 

7(0 

-i{C,fdv = -f{l3{C.)fdvo(3 
<;=^ dv = ^~‘^da. 


7(0 

7(/3-H0) 


dv 


The formula for P follows from (i)-(ii). 


□ 


The formula for P explains the need for the operator R in Lemma rra The 
proof of quantum ergodicity will be based on a reduction to the classical L^-mean 
ergodic for the symbol using Egorov’s theorem. In the case where (3 is ergodic, the 
mean ergodic theorem for T on Lf{B*Y, (,)) states that: 

hm ||l^T"(/)-P(/)||^0. 

N^ca Jy ^' 
n—1 

In the automorphism case, P^ay) = w(A) and we can deduce the desired quantum 
ergodicity result directly from the mean ergodic theorem. In the endomorphism 
case, P{aA) = ca;Neu(A) 7 , so we need to construct an operator R which contributes 
the factor of 7 . 

In the Dirichlet case, the roles of T and T* get interchanged. By essentially the 
same proof, we then have: 

Proposition 2.2. Consider the operator T* on Lp'{B*Y^da). We claim: 

• (i) The unique positive T*-invariant density is given by 7 “^. 

• (ii) T* is unitary relative to the inner product ((,)) on B*Y defined by the 
measure dv = j'^da. 
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• (iv) When (3 is ergodic, the orthogonal projection P onto the invariant L^- 
eigenvectors has the form 


P{f) = 


-1 


voi(s*y) 


[Jflda]-f ^ = ca;Dir(/) 7” 


B»Y 


3. Piecewise smooth manifolds 

Let ri be a bounded subdomain of M" with closure il. 

Definition 3.1. We say that C M" is a piecewise smooth manifold if the bound¬ 
ary Y = on is strongly Lipschitz, and can be written as a disjoint union 

r = u • • • u u s, 

where each Hi is an open, relatively compact subset of a smooth embedded hyper¬ 
surface Si, with Q lying locally on one side of Hi, and where Y is a closed set of 
{n — 1)-measure zero. 

The sets Hi are called boundary hypersurfaces offl. We call S the singular set, 
and write Y° = Y \ Y for the regular part of the boundary. 


Remark. There appears to be no standard definition of ‘piecewise smooth man¬ 
ifold’. Our definition is rather broad; we have chosen it since it seems to include 
all known examples of ergodic Euclidean domains. We remark that it includes, 
for example, nonconvex polygons, and arbitrary convex polyhedra which are not 
included in the class of manifolds with corners as considered in EEH], fZZ^ and 
m- However, the increase of generality is more apparent than real since the 
proofs in |ZZw| and |Bu| undoubtedly go through for piecewise smooth domains as 
defined here. Moreover, |ZZw| and |Bu| deal with Riemannian manifolds, not just 
Euclidean domains. 

The operators need to be properly defined on a piecewise smooth manifold. 
For all boundary conditions considered here, it is defined as the self-adjoint operator 
associated to a closed, semibounded quadratic form. The quadratic forms are 

(3.1) <3Neu(M) = [ \Vu{z)\'^dz ueH^{^), 

JQ 

for the Neumann boundary condition, 

(3.2) Qnoh{u) = [ |Vu(z)pdz-b [ K{y)\u{y)\'^ da{y) m e i7^(B), 

JQ JY 

for the Robin boundary condition, 

(3.3) Qx{u)= f \Vu{z)\‘^ dz-\- [ {Ku){y)u{y)da{y) u€H^[Vl), 

Jn JY 

for the d'^-Robin boundary condition, and 

(3.4) goi,.(M) = f \Vu{z)\'^ dz ueH^i^), 

Jn 

for the Dirichlet condition. In (lOl . K is a first-order classical pseudodifferential 
operator whose kernel is supported in Y° x Y°. Here closedness of (13.311 (hence 
self-adjointness of A^) requires that K be self-adjoint, while ellipticity of the 
boundary problem requires that K have non-negative principal symbol (the form 
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is positive when K is positive, and semi-bounded provided the principal symbol is 
non-negative). 

Suppose that n is a piecewise smooth manifold. We shall define a ‘compressed’ 
cotangent space T*Y = T*Y° U S which is the space of which our wavefront set, 
dehned below, will be a subset. For each boundary hypersurface Hi we define 
a compressed cotangent bundle T*Hi by identifying the fibre T*Hi above each 
y € dHi to a point; we give this space the quotient topology. Then we define the 
topology on T*Y by gluing together the T*Hi along their boundaries. The space 
T*Y essentially ‘ignores’ the frequencies over S. 

Let us now define certain classes of operators on L'^{Y), depending on a real 
parameter h £ (0, ho]. First, we define C°°{Y) to be the restriction of to 

y, and define C°°(Y'^) similarly. We define a residual operator to be one with a 
kernel Kh{y, y') (all kernels are written with respect to Riemannian measure on Y) 
which are smooth in {y,y') for sufficiently small h, and are 0{h°°) as an element 
of C°°{Y'^) (in the sense that each seminorm of Kh is 0{h°°)). We next dehne a 
semiclassical pseudodifferential operator of order m on T to be an operator which 
can be expressed as the sum of two terms, one a semiclassical pseudodifferential 
operator on y°, parametrized by h G [0, ho] for some ho > 0, whose kernel is 
properly supported in Y° xY°, and the second a residual operator. We shall work 
here only with the most classical types of operators, namely those which are either 
differential operators, or with symbols a{y, y, h) which are smooth functions of y, y 
and h, compactly supported in y. 

We similarly dehne a semiclassical Fourier integral operator (FIO) to be the 
sum of a semiclassical Fourier integral operator with kernel properly supported in 
Y° X Y°, and a residual operator. For a pseudodifferential operator or FIO Qh, the 
(semiclassical) wavefront set support WF'((5?i) is well dehned and is a closed subset 
of T*Y° X T*Y°; for pseudodifferential operators, it is a subset of the diagonal. 

We may dehne the operator wavefront set more generally. Let Th be some 
operator on Y, and let q,q' G T*Y. We say that {q,q') ^ WF'(r) under the 
following conditions: if q, q' G E, then if there are smooth functions (p, p' with 
(j){q) 7 ^ 0, (j){q') 7 ^ 0, and (j){y)Th4>{y') residual; if g G E, g' G T*Y°, then if there 
is (j) as above and a pseudodifferential operator Qh which is elliptic at q' with 
4>ThQh residual; similarly if g' G E, g G T*Y°; and for q,q' G T*Y°, if there 
are pseudodifferential operators Q,Q' elliptic at q,q' respectively, with QuThQ'h 
residual. The last condition is equivalent to the usual condition when Y is smooth. 
The operator wavefront set WF^(T) is then a closed subset of T*Y. 

The billiard map [3 : B*Y° 'f*Y is dehned on the open ball bundle B*Y° as 
follows: given {y,y) G T*Y, with \y\ < 1 we let (y, C) G be the unique inward¬ 
pointing unit covector at y which projects to (t/, y) under the map ^ T*Y. 

Then we follow the geodesic (straight line) determined by {y, () to the hrst place 
it intersects the boundary again; let y' GY denote this hrst intersection. If y' G E 
then we dehne P{y,y) = y'■ Otherwise, let y' be the projection of (f to T*,Y. Then 
we dehne 


Piy,v) = 

The map (3- : B*Y° j'*y is dehned similarly, following the backward billiard 
trajectory (that is, the straight line with initial condition [y, 2{f ■ i'y)vy — C))- 



QUANTUM ERGODICITY OF BOUNDARY VALUES OF EIGENFUNCTIONS 


13 


The billiard map is a symplectic, hence measure preserving, map with respect to 
the standard symplectic form on T*Y. This follows from the fact that the Euclidean 
distance function d{y,y') is locally a generating function for /3; that is, the graph 
of /3 in a neighbourhood of (yo, r]o,yo,Vo) i® given by 

(3.5) {{y,-"^yd{y, y'), y', ^y'd{y, y')) 

for {y,y') in a neighbourhood of (i/Oj2/o)- denote the graph of (3 by Cbuiiard, 

(3.6) ^billiard = graph/3 = {{I3{q),q) \ q S 7^^}. 

This is a smooth Lagrangian submanifold of B*Y° x B*Y°. For strictly convex fl 
it is given globally by for y,y' € Y°, but this is not true in general. This 

causes extra difficulties for nonconvex domains which are dealt with in Section [TTl 
It is clear that when S is nonempty, /3 does not map the open ball bundle to itself. 
We define 7Z^ C B*Y° to be /3-^(B*Y°) and, more generally, 7Z^+^ = for 

natural numbers k. Thus TZ^ consists of the points where /3* is well defined and 
maps to B*Y°. Similarly we dehne TZ~^ = f3Z^{B*Y°) and TZ~^~^ = /3Z^(7?.“^). 
Clearly TZ} Z TZ? , TZ~^ D TZ~'^ D ... and it is shown in |CFS| that each 7?.^ 
has full measure. Consequently, if we define TZ°° = HtTZ^, then /3 is a measure¬ 
preserving bijection on 7Z°°. 

4. Structure of the operators Eh and Fh- 

In this section, we shall decompose the operator Ft, defined by as well as 

the related operator Eh, into microlocal pieces. Here Eh is the operator with kernel 

Eh{y,y') = ‘2Go{y,y',h~^) 

y.y'eY 

where Go{y,y',X) is the free outgoing resolvent kernel (A — (A -b on R". 

The first two pieces shall be standard types of operators, namely Fourier integral 
operators, which are well understood. The third piece is a ‘left over piece’, but its 
operator wavefront set shall be quite localized. 

The reader is reminded that in the section we assume that H is convex. 

For simplicity, we first consider the case when Y is smooth. We denote the 
sphere bundle {{y,r]) G T*Y° \ \r]\ = 1} by S*Y°. 

Proposition 4.1. Assume that n is a smooth convex domain. Let U he any 
neighbourhood of S*Y° x S*Y°. Then there is a decomposition of Eh 

Eh = Eiy + E2,h + E^y, 

where Ei is a Fourier Integral operator of order —I associated with the canonical 
relation Cbiiuard given by (ESI), E 2 is a pseudodifferential operator of order — 1 
and £'3 has operator wavefront set contained in U. The principal symbol of E 2 in 
B*Y \ Ui (where Ui is the projection of U to the first factor) is i/ 7 , where 7 is 
defined by (CSl). 

Similarly, there is a decomposition of Fh as 

Fh = Fiy + F 2 .I 1 + Fsy, 

where Fi is a Fourier Integral operator of order zero associated with the canonical 
relation Cbiiiiard, F 2 is a pseudodifferential operator of order —1 and F 3 has operator 
wavefront set contained in U. 
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Proof. We first tackle the simpler operator Eh- The kernel of is given by 
C\^-^{\y - Ha«2_i(l2/ - y'l/h) 

The Hankel function Ha|^y 2 _i (0 is conormal at t = 0 and as t ^ oo, b{t) = 

expansion in inverse powers of t, with leading term ~ . 

We introduce cutoff functions 1 = (fiHy — y'\) + 4’2{\y — y'\/h^^^) + 4>3{\y — y'\,h), 
where (flit) is supported in t > to for some to > 0 to be chosen later, and (l) 2 it) is 
equal to 1 for t < 1 and equal to 0 for t > 2. (The power 3/4 in ^2 could be replaced 
by any other power strictly between 1/2 and 1.) Then, (j)i{\y — y'\)Eh{y,y') has a 
kernel of the form 

C/i-("-2)e4y-y'IA0^(|y _ y(\)bQy _ y(yh), 

where &(t) has an expansion in inverse powers of t as t ^ oo, with leading term 
~ This is manifestly a semiclassical FIO of order —1, and since fl 

is assumed convex, the phase function \y — y'\ parametrizes the billiard relation 

C^billiard- 

We next show that the kernel (f^Eh has operator wavefront set supported in the 
set U, if 6 is chosen sufficiently small. It is sufficient to show that if Qh € 4'°(y) 
satisfies WF'{Qh)oU = 0, then Qh^ifsEh) is residual, and similarly, if Qh G 
satisfies WF'((5;i) o [/ = 0, then {(f^Eh) o Qh is residual. The kernel of Qh o {(fsEh) 
is given by an oscillatory integral of the form 


h^- 


n 


E^y-y")-^/\iy'\ y, h)E\y”-y'\/^bi\y" - y'\/h)M\y" - y'\,h) dy" dy 


where <l{y",r]) is supported away from \ri\ = 1. The phase is stationary when 
y = dy"\y" — y\. Since (/a is supported in the region where \y — y'\ < <5, and 
\dy"{\y" — 2/1)1 ^ 1 as \y — y'\ 0, this means that the phase is never stationary 

if to is sufficiently small. Repeated integrations-by-parts show that the kernel is 
residual, since we gain an h each time we differentiate the phase and lose at most 
when we differentiate (f^. The computation for Qh is similar. 

Next we analyze the kernel (j) 2 Eh. We shall show that this kernel is pseudodif¬ 
ferential when microlocalized away from the set |? 7 | = 1. To do this, we write the 
kernel of (j) 2 Eh as the distributional limit, as e ^ 0, of the oscillatory integral 

(4.1) 2(2,)-"V->,(k^) 

We write ^ -I- fvy', where S TyiY, ^ G R and Vy' is is the inward pointing 

unit normal at y'. The kernel can be written 


2(27r)-"/i2-"02( 


\y-y'\ 

/l3/4 


gdy-vf-ffh 


gi{v-v')-t'y'Clh 


dC" dC- 


To localize away from \y\ = 1, we introduce cutoffs 1 = 

where ifiit) is supported in t < 1 —ti, ' 4 ’ 2 .it) is supported in 1 —2ti <t< l-|-2ti and 
ipzit) is supported in t > 1 -I- ti. Inserting the cutoff ipi means that — 1 < 0, 
so we can perform the f integral using the formula 


1 

27r 




e - {a + z0)2 


dC = 


ze 


ika 


2a 


k > O.a > 0 



QUANTUM ERGODICITY OF BOUNDARY VALUES OF EIGENFUNCTIONS 


15 


to get 


(4.2) z(27r)-"+i/i2-"02( 


\y-y'\ 

/l3/4 


Av-v )•«" /h £- -(^ 11 ) d^W . 

v'l-IC'ip 


As h 0, {y — y') ■ Vy/ jh ^ 0 on the support of ^ 2 , since \y — y'\ < 2 ft, 3/4 on the 
support oi (f >2 but {y — y')-Vy = 0{\y — y'\'^). Hence we may expand the exponential 

e* 


od!/ V )-‘^y'V^ Taylor series centred at zero: 

w-i _ 

(4.3) |e*(J/-^')-vVi-l?''IV/*_ ^ -(*(y-y')-V\/l-|?"PA)^ 


3=0 


<Cn 


{y-y') 


N 


Consider one of the terms 

(4.4) ^(2^)-"+1 /i2 - J e^(v-v')<"/h(^^^y _ y/) . - |^ll|V/i)Vi(ell) 


Since {y — y') ■ Vyi = 0{\y — y'\^), after 2j integrations by parts we eliminate the 
vanishing at the diagonal and gain a factor of . This is therefore a pseudodiffer¬ 
ential operator of order — j. Multiplication by the (j )2 factor only changes this by a 
residual kernel, since H4.4II vanishes rapidly as\y — y'\/h ^ oo. Similarly, the error 
term is becoming more and more regular. Therefore 63 is a pseudodifferential 
operator of order —1. Moreover, we see from the form of that the principal 
symbol of this operator is i/— jiyp = i/l- 

Similar reasoning applies to the cutoff ^ 3 , using instead 


1 

2 ?! 


p»^C. 


C + a‘ 


■dc 


g-ka 

——, k,a>Q. 
2 a 


This gives the kernel 


(271) 


— n +1 T2— n 






-(y-y')'V 


vW^ 


7/’i(C")d^ll. 


We can similarly expand the exponential in a Taylor series to show that we get a 
pseudodifferential operator of order —1. 

The operator 63 with cutoff '02 inserted has operator wavefront set arbitrarily 
close to S*Y° X S*Y°, and hence within U provided ti is sufficiently small. This is 
shown as for (j) 2 Eh, with the help of Theorem 7.7.1 of |Hd| which gives e-independent 
estimates on all seminorms of the composition of the operator with Q (on the left) 
or Q (on the right). Hence, if we define Ei to be the operator (j)iE, E 2 to be 
the operator with cutoffs 02 (' 0 i + V'a) and E 3 to be the remainder, we have a 
decomposition which satisfies the conditions of the theorem. 


To deal with the operator Fh, we argue similarly. Using the same (pi cutoffs as 
before, the operator piFh is an FIO of order 0 and the operator p^Fh has wavefront 
set contained in U. To deal with the remaining term, we write the kernel of p 2 Eh 
as the distributional limit, as e —> 0, of 


2(27r)-"/ii-" 


d^iv-v'X/h 


, f\y-y'\ 

d3/4 



(4.5) 
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We decompose ^ + (^Vyi as before, and write the kernel 

(4.6) 




-< 




C^ + IC'IP-I-' 

We use the cutoffs 1 = V'i(C'') + V' 2 (C'') + as above. Inserting the cutoff V'l 

means that — 1 > 0 , so we can perform the C (oscillatory) integral using 

\ /■ ^ika 




27ri 


c 

A2 / i ■n\2 O ’ A), Ci > 0 

C — (a + zO)"^ 2 


to get 




Following the reasoning above, this appears to be a pseudodifferential operator with 
symbol 1. Similarly, the term with ips appears to be a pseudo with symbol 1, which 
would give us the identity operator modulo an operator of order —1. However, 
the identity term is not present in the kernel Fh since the kernel of the identity is 
supported at the diagonal and does not appear in the restriction of the kernel of 
du^,Go{y,y') to the boundary. In fact, it is the ‘same’ identity operator that turns 
up in the jump formula for the double layer potential in Thus, this piece of 

Fh turns out to be of order —1, as for Eh- 

Finally, as for Eh, the operator with cutoff rp 2 may be shown to have wavefront 
set arbitrarily close to As*y, and hence within 17 by a suitable choice of the cutoffs 
ipi- Hence, if we define Fi to be the operator f/fiF, E 2 to be the operator with 
cutoffs 02(V’i + V’a) and F 3 to be the remainder, we have a decomposition which 
satisfies the conditions of the theorem. □ 


Remark. A related analysis of Fh is given in m (it is denoted N{k + ir) there). 

A more complicated version of this Proposition is valid when dfl has singularities. 
We first define 

(4.7) E = S*Y°Ui:cf*Y, 

and recall the notation TZ^ from Section |21 We denote the complement of a set S 
by 5'. 

Proposition 4.2. Let Ll be a convex, piecewise smooth domain. Let U be any 
neighbourhood of'E.x {Ti})' U {TZ~^y x S. Then there is a decomposition of Fh 

Fh = Fi^h + F2^h + F^^h, 

where Fi is a Fourier Integral operator of order zero associated with the canonical 
relation Cbiiiiard, F 2 is a pseudodifferential operator of order —1 and F 3 has operator 
wavefront set contained in U. 

Proof. We choose a function (j) on Y so that <() = 1 on a neighbourhood of E, so 
that 

{iq,q) I 7r((7),7r(g') G suppef} C U, 

{(/3(g), 9 ) I 7r(g) G supp(/)} C U, 

{(g,/3"^(g)) I 7r(g) G supp(/)} C U. 
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We may write 

F = + (t>F{l -cj)) + {l- (t>)F(j) + (1 - (f)F{l - (j)). 

The term (1 — (j))F{l — (j)) is supported away from the singular set in both variables, 
so may be treated by the argument above. We claim that the remaining terms have 
wavefront set contained in U. This is clear for (jiFcj), and the argument for the other 
two terms is similar, so we concentrate just on (1 — We will show that the 

wavefront set is contained in 

(4.8) {iq,q') I TTiq), 7 T{q') e suppc/)} U {(/3(g),g) | Tr{q) G supp(?!)} 

Thus, let {q, q') be a point not contained in 114.811 . We want to show that the kernel 
(1 — (j))F(j) is regular at {q,q'). If T^{q') ^ suppt)), this is obvious, so assume that 
7 r(q') G suppi^, 7 r(g') ^ supp()), and that q = (yo^Vo) ^ PW)- Then, there is a 
smooth function a{y,r]) with 7 ( 5 ) = 1 , a smooth function (p{y) with (p(Tr{q')) = 1 , 
and with the support of a disjoint from 

{P{q') I 7 r(g') G suppc^}. 

Let g 1 be a semiclassical pseudodifferential operator with symbol a. Then the 
composition A{1 — (j))F(f>(j) is represented, modulo a residual term, by an integral of 
the form 

(4.9) j F^y-y''^-'^/^F\y"-P\/'^a{y'\y)r{y'',y){(j)p){y')dydy''. 

Here r is a smooth function, since on the support of ap, y" and y' are separated. 
The phase in this integral is never stationary on the support of the symbol, by 
construction. Hence this operator is residual, proving that {q,q') ^ WF'((1 — 

<t>)FP). □ 

In order to deal with products of operators involving T 3 or F^ we need the 
following proposition. 

Proposition 4.3. Assume that the semiclassical operator Ah is either a pseudodif¬ 
ferential operator or a Fourier integral operator with compact operator wavefront set. 
//WF'(H)oWF'(F 3 ) = 0, then AF^ is residual. Similarly, z/WF'(F '3 )oWF'(H) = 0, 
then FfA is residual. 

Proof. We first observe that if Rh is a residual operator, then WF'(i?F 3 ) is con¬ 
tained in {{y,y) \ y G S}. Indeed, this property for F follows readily from the form 
(tra of the operator, and it holds for Fi and F 2 since these are FIOs. Therefore it 
holds also for F^ = F — Fi — F 2 . 

To prove that AF^ is residual, notice that the sets 

D, = {q'Gf*Y\ 3{q,q')GWF'iA)}, 

D 2 = {q'Gf*Y\ 3 (g',g")GWF'(F 3 )}. 

are closed, and they are disjoint by hypothesis. Hence there exist disjoint open 
sets Oi D Di, O 2 A D 2 . Choose Qh G such that the symbol of Q is I in 

Oi X (0, ho) and zero in O 2 x (0, ho) for some hg > 0, and so that the kernel of Q is 
supported away from S in both variables. Then QF 3 is residual by hypothesis and 
and H(Id—Q) = A — AQ is residual by the symbol calculus for pseudodifferential 
operators, and supported away from S in both variables. Writing 

AFs = AiQFs) + (A(Id-Q))F 3 , 
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and using the observation above we see that AF'^ is residual. The result for F^A 
follows by taking adjoints. □ 

We conclude this section with a crude operator bound on the operator Fh- 

Proposition 4.4. Let ft he a convex domain with corners. Then the operator 
norm of Fh admits a bound 

( 4 . 11 ) \\Fh\\mY)^mY)<Ch-(^-^\ 

Proof. The kernel of Fh{z, z') is given by 

h~‘^Vy ■ (z - z')f{\z - z'\/h), 

where f{t) is symbolic as f ^ 0, with f{t) = + g{t), where g{t) = 

as f —> 0, and bounded as t ^ oo. Thus, the kernel is given by 

■ |/_ - z'\lh). 

The first kernel is bounded on Lf by the theory of singular integrals on Lipschitz 
submanifolds; see IlmH!. The Lf operator norm of the second kernel, whose sin¬ 
gularity at the diagonal is 0{\z — may be crudely bounded by Schur’s 

Lemma, giving a bound □ 

We remark that under any decomposition as in Proposition 14.21 Fi and F 2 are 
uniformly bounded on L^, so the bound in is also valid for F 3 . 

5. Local Weyl law 

In this section we shall prove Lemma o for the Neumann boundary condition. 
Let us first prove (truii for Ah = Id, which is the statement that 

2vol(r) 


(5.1) 




E 

A,<A 




rol(n) 


We use the Karamata Tauberian theorem, which states (see eg na, p89) 

Proposition 5.1. (Karamata Tauberian Theorem) If g, is a positive measure on 
[0,oo), and a> 0, then 


e *^dg{X) ~ at 


t^O 


implies 


/o r(Q! +1) 


To prove eu we consider 




e{t) = E' 

1=1 

This is equal to the trace of the operator E{t) from (112.1011 in the appendix. By 
dma, we see that 

e{t) ^ 2(47rt)“"’/^ vol(y) as t > 0, 
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SO by Proposition 15.II 




2 voi(r) 


Xj <CA 


On the other hand, 


NiX) 


(47r)"/2r(n/2 + 1) 
vol(O) 


A". 


-A" 


(47r)"/2r(n/2 + 1)' 

so follows. 

Next we prove the lemma for zeroth order pseudos Ah supported away from E 
in both variables. We proceed through a series of reductions. 

(i) First, if Ah is in (that is. Ah = hAh for Ah € 4'*’°(F)), then 

( 5 . 2 ) ^ 1 ^^^ 


Xi<x 


This follows easily from and the fact that the operator norm of Ah is uniformly 
bounded in h ( |DSj| , Theorem 7.11). 

(ii) Second, if the support of the symbol of A at /i = 0 is contained in {|^| > 1}, 
then we also have (EJ. 

To see this, we use the fact that Fh^u'^ = uj, to write the left hand side of (15.211 


as 


lim 
A —*00 


1 

Mx) 


E {F*h^Ah^Fh,u],u^^. 

Xj<X 


We use Proposition 01 from the previous section to analyze F^AhFh- Due to the 
condition on the symbol of A, we can find an open set U as in ProDOsition l4.2l disioint 
from WF'(A). Let F = Fi + F 2 + F^ he a decomposition as in Proposition l4.2l with 
respect to U. Then WF'(A)oWF'(Fi) = WF'(A)oWF'(F 3 ) = WF'(FC)oWF'(A) = 
WF'(F'I) o WF'(yl) = 0. Thus by Proposition 14.HI the operators 


A{Fi + Fs) and (Fi + FsTA 

are residual. Hence, using also Proposition l4.4l we have an 0{h°°) estimate on the 
operator norm of F* AFj unless i = j = 2, which certainly gives us the required esti¬ 
mate in these cases. In the remaining case, i = j = 2, F 2 AF 2 is a pseudodifferential 
operator of order —2, so the required estimate follows from (i). 

(iii) Third, there exists an integer k such that, if the symbol of Ah at /i = 0 is 
sufficiently small in norm, then 

(5.3) lim sup -^1 E \ < e. 


To prove this, note that by (i) one may assume that the symbol is independent 
of h and by (ii), that the support of the symbol is in {|^| < 2}. Then, the norm 
of the symbol, for sufficiently large fc, controls the operator norm of Ah, uniformly 
in h. Thus, if the C* norm is sufficiently small, the right hand side is smaller than 
e. 

(iv) Since the symbols of differential operators are dense in the topology in 
the space of symbols supported in {|^| < 2}, it is enough, by (iii), to prove the local 
Weyl law for differential operators. 
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(v) Next we note that for odd order monomial differential operators A, as above, 
the limit is zero. To see this, note that 


1 


I] = 

\,<\ 


1 


x,<\ 


but A + A* is in d'*’ ^(T) if A is an odd-order monomial, so by (i), the limit is 
zero. This clearly agrees with so it remains to treat even-order operators. 

(vi) Finally, every even-order differential operator is the difference of two positive 
differential operators of the same order. Hence we may restrict attention to positive 
even-order operators supported in a single coordinate patch. 


Thus, let Ah = where P is a positive differential operator of order 2k. 

Consider the quantity 

OO 

eA{t)=Y,e-^^hf^{Pu),u)). 

j=i 

Let us consider which is given by 


^^'i{Pu],u]) =tTPE{t). 


By (I12.1dll . this is given by 
(5.4) tvPE{t)= 2 


(47r)"/2 vol(S'" 


J a(P)d/rNeu + 0(i-^-"+^). 


B*Y 


Integrating k times in we find 

eA{t) = ^ 0 . 

B-Y 

Then 11.1511 follows from Proposition lb.ll as before. This completes the proof of the 
Lemma. 


6. Egorov Theorem 

In this section we prove a generalization of Lemma ll.41 First, we need to compute 
the principal symbol of the FIO Pi from Proposition ^21 

Proposition 6.1. Let Pi and U he as in Proposition 14-^ Then the principal 
symbol of Fi at {P{q),q) G Cbuiiard \U is 

a(Pi(/3(g),g)) 

where t is an eighth root of unity and \dq\ represents the symplectic density on 
B*Y. 

Remark. The Maslov bundle over Cbmiard is canonically trivial, since Cbmiard is 
the graph of a canonical transformation, so we shall ignore it. 
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Proof. We begin with the explicit expression (II.till for the free resolvent on R". At 
the billiard Lagrangian, y and y' are distinct, so we need the asymptotics of the 
Hankel function as it argument tends to infinity. This is 

OO 

i=o 

where oq = ■\/ 2 / 7 r. Moreover, the expansion can be differentiated term by term. 
Thus, the principal symbol of F{X) at the billiard Lagrangian (up to an eighth root 
of unity) is the same as that of 

(27^)-(”-l)/2A("+l)/2g^A|y-y'| |\v - y'\. 

We see from this expression that the operator has semiclassical order zero, and its 
symbol is |Ho4j . |PUrj 

\y - \y - y'\\dydy'\^^'^. 

It is more geometric to give the symbol in terms of the symplectic half-density 
\dydr]\^^'^ = \dy'dri'\^^'^ = \dq\^/'^ on T*Y. To do this, we express dy' in terms of dry, 
keeping y fixed. Since rji = dy^ ly — ry^|, we find that 

|dry| = det \ {y + sief) - (y' -b tie')| j |dy'|, 

where is an orthonormal basis for TyY, and e) an orthonormal basis for Ty/Y. To 
make things a bit clearer we first compute this in the two dimensional case, n = 2 . 
We choose coordinates so that y is at the origin, y' = (0,r), ei = (cos a, sin a) and 
62 = (cos/3, sin/3). Thus, we are trying to compute 

dsdt 

This is equal to 

cosacos/3 = |y - y'r^a,.Jy - y'\d^^,\y - y'|, 
so this gives the symbol 


I (s cos a — t cos /3, s sin a — r — t sin fd) |. 


( 6 . 1 ) {d„^,\y-y'\)\y-y'\^^^{\y-y'\ ^d^Jy - y'\d^^,\y - y'\^ 


^^^Idydryl^^^ 


1/2 1/2 _ / i [ y ', v ') ^l/^ 

Adzzjy-y'l/ \ y d\ \^(f3(y' ,7]'))) 


dydr] 


1/2 


where a is defined in C3- 

In the higher dimensional case we compute Idryj as follows: let X C M" be the 
subspace of R” 

x = {TyY)niTy,Y)nij^y„ 

where ly^y' is the line joining y and y'. We will initially assume that all intersections 
are transverse, so that X has codimension three. Then TyY n X has codimension 
two in TyY, and similarly for TyiY . We choose Cg = 63 ,..., e(j_i = e„_i to be an 
orthonormal basis of X. We choose ei G TyY to be in the plane of ly^y/ and Uy, 
and 62 to be orthogonal to both ei and ly,y', and choose e'l and 62 similarly. Then, 
61 ,..., e„_i is an orthonormal basis of TyY, and e[,..., e'„_i is an orthonormal 
basis of Ty^Y. 
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Let US choose coordinates in R" so that X = {zi = Z 2 = Z 3 = 0}, with y = 
(0,..., 0), y' = (0,0, r, 0,..., 0). Without loss of generality we may assume 

ei = (cos a, 0 , sin a), 

e2 = (0,1,0), 

e\ = (cos 7 cos/3, cos 7 sin/3, sin 7 ), 
e '2 = (—sin/3, cos/3,0), 


where we write only the first three coordinates since the rest are zeroes, and ej = e) 
is the j + 1st unit vector in R". Thus we need to compute 


( 6 . 2 ) 


det 


ds’-dP 


(/(si,. 


/ = |(/i,---,/n)|, and 


,Sn-i,ti,-.-,tn-i)), where 
fi = Si cos a — ti cos 7 cos /3 + t 2 sin /3, 
/2 = S 2 — ti cos 7 sin /3 — t 2 cos /3, 

/a = Si sin a — r — p sin 7 , 
fi = Sj-i-tj-i, j > 4. 


A simple but tedious computation shows that ijOl is equal to r cos a cos 7 . If 
the transversality assumptions are not true, then this result follows by perturbing 
to a case where they are satisfied, and taking a limit. Thus, in general the symbol 
is given by EIJ. □ 


Remark. The proof also shows that the principal symbol of Ei at (/3(g), g) G 
C^biiiiard \ U IS equal to 


(6.3) 


<^{Eh){P{q),q) 


—ir 

\/i{l3{q))i{.q) 


Mgr/^ 


Lemma 6.2. Suppose Ah is a pseudodifferential operator of order zero, withWF'(A) C 
TZ~^. Then, one ean express 

{F*^fAhFl^ = Bh + Sh, 


where ||>S'ft,||i 2 ^i ^2 < Ch, and Bh is a pseudodifferential operator of order zero with 
symbol 


(6.4) 


a{B) 


7 (g) 7 (/ 3 ^(g)) g G 

0 , g ^ 7^^ 


Proof. Suppose that a neighbourhood V of (S x S)U(S x (7?.“^)') U (7^' x E)U As*y 
is given (we shall choose it later in the proof). We decompose F = Fi + F 2 + F 3 as 
in Proposition 01 with respect to V. Consider the product 

(6.5) (7^1* + Ff + FffA{F, + F 2 + F^)^. 

Consider terms in the expansion which contain at least one F 3 or Ff (it will 
suffice to consider just F 3 ). These contain a string of the form 

( 6 . 6 ) ... A{Fi or F 2 )(Fi or F 2 ).. .{Fi or ^ 2)^3 ... 

The substring A{Fi or F 2 ){Fi or F 2 )... (Fi or F 2 ) = G is an FIO associated to 
the canonical transformation (3^, for some k with 0 < k < N. Thus, its operator 
wavefront set is contained in 

(6.7) Uo<fc<^{(g,/3-'=(g)) | g G WF'(A)}. 
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By our assumption on the wavefront set of we have (3 ^{q) well defined for 
q G WF'(A). We now choose F to be a neighbourhood of 

(6.8) S X (7^^)'U X S 

such that no points of the form {q, *) are in V, for q G WF^(A). This is possible 
since (lira is compact, and the set WF'(A) x T*Y is a closed set which is disjoint 
from (16.811 . Then the FIO G is such that WF'(G') o WF'(F 3 ) is empty, so by 
ProDOsition Id.dl the composition (16.611 is a residual operator, and hence is bounded 
on with an 0{h°°) bound. Since, by Pronosition 14.41 we have a bound 
on all the remaining terms in the product (16.611 . it follows that each of these terms 
has an 0{h°°) bound on its operator norm, hence certainly are of the form Sh 
above. 

Thus, it remains to consider factors which do not contain an or F^. Consider 
next all terms in the expansion of (ES which contain only Fi ’s and F 2 ’s (including 
their adjoints), and at least one F 2 . Since FIOs of order zero associated to a 
canonical transformation are bounded on L^, and F 2 is order —1, all these terms 
are bounded on with an 0{h) bound, so these terms are of the form Sh above. 

We are thus reduced to analyzing the term (F^)^BF(^. We have shown that 
Fi is an FIO whose canonical relation is Cbmiard = gi'aph(/3) with symbol given 
by Proposition 16.11 and hence, its adjoint F^ is an FIO with canonical relation 
graphand with symbol 

= CT(Fi)(/3(g),g) ^ {/3{q),q) i V. 

Thus, by the calculus of FIOs, the composition {F*)BFi is a pseudodifferential 
operator with symbol j{Piq))~^jiq)b{l3{q)). By induction, we see that the term 
{F*)^BF-f^ is a pseudodifferential operator with symbol (16.411 . 

This completes the proof. 

□ 


7. Proof of the main theorem — Neumann boundary condition 

In this section we prove ca for the Neumann boundary condition. It is sufficient 
to show that for any self-adjoint A, 

(7.1) ((AuJ,?!^) - WNeu(^))^ = 0, 


(This gives us the limit (11.511 for a density one subsequence depending on A; a 
standard diagonal argument then shows that this dependence can be removed.) As 
explained in m there are extra complications due to the fact that n~TTll is not 
an automorphism. The main novelty in the proof beyond consists in how to get 
around them. 

The proof consists of several steps. First we show that 


(7.2) 


lim 
A— KX> 


1 

mX) 


Y {iM.iD 

Xj <A 


^Neu(A) .. 
u;Neu(Id)" 


2 

= 0 . 
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(Recall that is the normalized boundary trace: = uJ/||M^||i 2 (y).) It is suffi¬ 

cient to show that for any e > 0 , we have 


(7.3) 


lim sup 

A—»oo 


1 

mi 


E ((^“S.«S> 

Aj <A 


‘^Neu(-^) 

‘^Neu(Iti) 


iI'w5iil2(y)^ 


< Ce, 


where C depends only on fl and A. Notice that in GS, exactly one of the uj’s is 
L^-normalized. This L^-normalization is a nuisance, but is not surprising, because 
in our proof of we make no reference to interior values of the eigenfunction, 
and the norm of in is dictated by the size of the eigenfunction in the interior. 
In the second step, we show that 

(7.4) ||M^|||2(y) —> WNeu(Id) 

along a density one subsequence. Then G3 and (El together imply the result of 
Theorem 11.11 


Step 1: Proof of 117.311 . We first choose a self-adjoint pseudodifferential operator 
R whose principal symbol r is a good approximation to cy, in the sense that 


(7.5) 


7 <e, 


vol(S'"-i)vol(II) 

4vol(B*y) 


We choose a real symbol a" such that a" vanishes at the boundary of B*Y^ and we 
have 

^Neu(c: ) — ^Neu(^)5 ^Neu((c: CL ) ) ^ 6. 

(Here we abuse notation by writing WNeu(a) for WNeu (^)5 etc, but this should not 
cause any confusion.) We then define 


(7.6) 


AN) 


N 








an averaged version of the symbol in El- By the mean ergodic theorem, we have 


7 ^{q)a"^^\q) ^ — w. [ 7 = cwNeu(a") = cwNeu(24) in L^. 

VOTiS Y ) Jb*Y 

(Note that we cannot make this claim for a since 07 “^ may not be L^, due to 
the simple vanishing of 7 ^ at the boundary of B*Y\ this is the sole reason for 
introducing a". ) Here c is as in El- We choose an integer N such that 

(7-7) || 7 fo^( 7 "^a"^^^ -ca;Neu(24))||^2(^.^j < e. 

We then choose a self-adjoint pseudodifferential operator A' whose symbol a' is 
close to a, and with additional properties required by Lemma 16.21 Namely, we 
require that WF'(H') is contained in n TZ~^, and disjoint from As*y, and is 
such that 

(7.8) (cUNeu)^!) CUNeu)^! )) < C, UJNeu((a U ) ) < C. 

Define analogously to in (17.611 . Then writing 

7_1/2(^,(^) -cyWNeu)^)) =7“fo^(a'^^^ - -C7WNeu(^)), 


and applying the elementary inequality for real numbers 
(7.9) (oi -(-----(- On) < -(-----(- a^) 
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with n = 2, we find 

(7.10) - C7WNeu(^)) II 2 < 4e. 


We then investigate the quantity 


(7.11) 


X! \{i^- ^Neu{A)R)u^j,U^j) 

A —*-00 (.-Aj 


A,- <A 


To work with this, we want to replace A with A! since A has been constructed so 
that Lemma 101 can be applied. Using 17.911 with n = 3, we bound H7.11|l by 


(7.12) 


Slim sup TTTTT 

A{\) 


E |((^-^') 




< E |((^' “‘^Neu(^')-R)w?.w5) 

I Aj<A 

+ E |((wNeu(^) - WNeu(^')).R'“j>'“' 


\j<C.\ A j A 

In the first term of (tTdUl . we may replace A'^ by 


N 




k=l 


since the u’’ are invariant under Fh^ ■ Also using the inequality 


(T<,*) <(T^u,u) 


for self-adjoint operators T, we can bound 17.1211 by 


(7.13) 


E - WNeu(A')i?)^u 5 ,u 5 ) 

A^oo ^ 


+ E ((^ - «?) + (‘*^Neu(A) - U;Neu(A'))' E “ 


A,- < A 


A1 < A 


We now use the Local Weyl Law, that is, Lemma to estimate the three lim 
sups in (Tmii . The second is bounded by e by rm . The third is given by 


4 (u;Neu(A) - u;Neu(A'))^ r r{qY 
vol(S'"-i)vol(U) Jb-y liq) 

Using arm and with n = 2 we can bound this by 


r (r(g) - C7(g))^ (C7(g))^ 

Jb^y 7(g) 7(g) 


where C depends only on Q. Using we see that this is bounded by Ce^ -|- C't. 

To estimate the first term, we use Lemma lti.2l and express + 

S, where S has an 0{h) operator norm bound and is a pseudodifferential 
operator with symbol given by a'^^\ By (EH) we may neglect the operator S, and 
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replace by Now, using the local Weyl Law for the operator — 

a;Neu(A')i?)^, the first term is given by 


(7.14) 


a;Neu(AW -CCNeu(A')i?)" =c / ■ 

J B*Y 


/W 


- WNeu(A')r)" 


dq 


< 2C 


[ Jb-y 7 Jb^y 7 



The first term is the key term which is controlled by (mnii thanks to classical 
ergodicity of the billiard flow. The second term is controlled by 113. This proves 
that is bounded by a constant times e. 


Now fra follows by replacing A in fTTTll with a;Neu(A).Id, and then subtracting 
the two expressions: 

r 1 \lf A WNeu(A) , 2 

hm sup -; 777 yy > ((A- 


A—^oo 


N{\) 


\j<\ 


‘^Neu(Id) 


(7.15) 


<21im sup ^ I {{A-uj^^^{A)R)u],u]) 

A^oo 


(( U;Neu(A) _ 


^Neu(Id) 


which we have shown is bounded by Ce, where C depends only on and A. This 
completes the proof of 113. 


Step 2: Proof of 17.411 . To prove this, we need to relate the boundary values of 
the Uj to their interior values, and use the L^-normalization of the uj on fl. We do 
this via Green’s formula: 

(7.16) Uj(z) = J^Gh^{z,y)u''j{y)da{y). 

Here, Ghiz,y) is the normal derivative, in y, of the free Green function for (—A — 
{h~^ + i0)^)“^ on K". (Thus, Fh is simply the restriction of this kernel to the 
boundary in z.) Let ip he a smooth nonnegative function compactly supported in 
the interior of hi. Then 

(7.17) = {Gl^'4f^Gh,u],u^^Y- 
We need to analyze the operator G'pip'^Gh- 


Lemma 7 . 1 . Let (p be a smooth function onY = which is identically one near 
S. Then 

(i) — (p) is a semiclassical pseudodifferential operator of order zero on 

Y; and 

(a) For each 1 < i < m, Gf^fPGh is bounded on LffFli), uniformly in h. 

(Hi) The limit o/a;Neu(G*V'^G(l — (p)), as the support of (p shrinks to Y, is 


1 


ip^{z) dz. 


(7.18) 


vol(fl) 
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Proof, (i) The kernel of the composition {tj}Gh)*{4’Gh){^ ~ 4‘): a function of 

( 2/1 y') X Y, is given by an integral of the form 

(7.19) 

;j-(n-l)( 2 ^)-(n-l) f h)f:^{z){l - <l)){y') dz. 

Jn 

Because y, y' are restricted to Y, and z is restricted to the support of if which is 
in the interior of fl, we do not see the singularity in the functions \y — z\, \y' — z\, 
or in a. The derivative of the phase function $(?/, y', z) = \y — z\ — \y' — z| in z is 
nonvanishing unless y = y'. Due to the cutoff function 1 — (j){y')i which is supported 
in the regular part Y° of the boundary, the kernel is smooth and order unless 
y = y' and both lie in a compact subset of Y°. 

To analyze further we write ^{y,y',z) = “ y')iVi, where r]i{y,y',z) is a 

smooth function, which we may do since $ vanishes whenever y = y'. Moreover, we 
have rji{y,y,z) = dy,^{y,y', z)\y>=y = I\.y{y - z/\y - z\), where liy is the projection 
onto the tangent space at y. Letting s = \y — zj, we find that 

drj = 7(7/, rf)dz, when y = y'. 

Hence we may change variable from z to (s, rf) in a neighbourhood of the set y = y', 
y, y' G supp 1 — (/), z G suppi/;. There, the kernel takes the form 

j e^^y-y'^-^biy, y', r,, h) drj, where & is 

CXD 

/ s^ ^ ds 

a{y,z,0),d{y,z,0)'tlj'^{z){l - (j){y')) , z = z{y,s,r]), 

0 

which is the kernel of a semiclassical pseudodifferential operator of order zero. 
(In the equation above, z = y + s9, where 9 G 5"“^, 9 ■ Vy = r], 9 is inward pointing 
at y.) Moreover, we see that the symbol b is nonnegative, and has the property 
that 

(7.20) WNeu(B/j) increases as either if oi \ — (f increases. 

(ii) This follows from (i). Let iJi,..., Hm be an enumeration of the boundary 
hypersurfaces of D. Each Hi may be embedded in an open smooth submanifold 
Hi. Let X be a function on Hi with compact support, which is identically 1 on 
Hi. Then, by (i), G'ftp'^GhX is a pseudodifferential operator of order zero, so is 
uniformly bounded. This implies a uniform bound M on ipGh acting from L'^(Hi) 
to Lf{Vt). Then mfM is a uniform bound for G^ijP'GhX on Lf{Y). 

(iii) To prove this, we analyze the integral ir7Tl?ll further. We are only interested 
in the principal symbol of G*'0^G(1 — (/>), so we may substitute the value of |a| at 
h = 0, which is 

Hy, z, 0)1 = - z|-i”-i)/ 27 (y, z) since j{y, z) = -p—^ • i^y. 

2 \y-z\ 

Thus, the principal symbol is given by 

00 

= i J s"(""^i 7 (y,z)V'^(z)(l - (()(y))s’""^(is, z = z{y,s,'n). 

0 
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Thus, WNeu(G*'!/'^G'(l — cf))) is given by 


1 

vol(r2) 


oo 

J J (j){y))-^{y,z)ds 

B*Y 0 


dydr] 

i{y,v)' 


The factors of 7 cancel. Now we change variables of integration to z and 9 G 5'"“^, 
where z = y + s9 and rj = 0 ■ Vy. Then, drj = -ydO, and for a fixed 0, dz = 'jdyds, 
so dy ds dr] = dz dO and the integral becomes 

voi(s-‘)voi(n) I ">***■ 

Clearly as the support of 4> shrinks to the singular set, this integral converges to 

(mini . □ 


We shall use the relation (m . Lemma o and the following result to deduce 

(TTll . 

Lemma 7 . 2 . Let {an) be a seqence of complex numbers. Suppose that for every 
(5 > 0 , there is a set of integers Ss of density at least 1 — (5 such that, along it, the 
oscillation of the sequence {an)neSs *5 most S. Then there is a set T of density 
1 such that the sequence {an)neT converges. 

Remark. The oscillation of a sequence ( 6 „) is defined to be 

lim sup \hni-hn\. 

m,n>N 


Proof. By hypothesis, for each n there is a set of integers 5'„ of density at least 
1 — 2 “" such that the oscillation of the corresponding subsequence is at most 2 “". 
By replacing Sn by U {1, 2,..., K} for suitable K, we can ensure that 


#gn{i,2,...M} 

M 


> 1 _ 2-"+i 


for all M. 


Consequently, Tn = S'„(n5'„+i(n... has density at least 1 —2“"’+^, and the sequence 
(afc)feGT„ has zero oscillation, that is, is a Cauchy sequence, hence converges to some 
number An. Since Tn (H Tm has positive density, for m > n > 3, we conclude that 
the An are all equal to some fixed A. 

By replacing by T„\ {1, 2,... K} for suitable K, we can ensure that |afe — A| < 
2“" for all k GTn. Assuming this condition, then, the set T = UnTn is density one, 
and (an)neT converges. To prove this, let n be given. Then choose Li so that 


sup |a/c — A| < 2“" for k G Ti, k > Li, 1 < i < n. 
Let L = maxi<„ Lp, then 

sup \ak — A| < 2“" ioY k G V, k > L, 1 < i < n. 


This proves that the sequence (a„) converges to A along T. 


□ 


Thus, by the Lemma, it suffices to show for each ^ > 0 that there is a sequence 
Ss of density at least 1 — CS such that the oscillation in the sequence 
along Ss is at most C6, for some constant C depending only on Q. 

Returning to im . we write = S and consider the equality 

(7.21) = ||Mj||i2(y.)(5'(l - (l>)u],ufiY + {S(j)u],u])Y. 
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By Lemma rm the operator S{1 — (j)) is pseudodifferential of order zero. Thus, by 
(ESI), we see that (S'(l —has a limit, namely vol(0)(vol(y)) ^WNeu(<5'(l — 
0)) > 0, along a density one subsequence. Now we choose ^ so that 0 < < 1, 

and ■)/) = 1 on a set of measure at least (1 — vol(r2). Then, by the local Weyl law 
in the interior of 

^ ^ a, <a 

On the other hand, each individual term {ip'^Uj,Uj)Q lies between 0 and 1, so we 
conclude that there is a subsequence Sg of density at least 1 — d on which the 
oscillation of the sequence is at most 25. 

Having chosen ip, we now choose (p with support so small that 


where tti?M is as in the proof of (ii) of Lemma o This is possible by and 

the Karamata Tauberian theorem. It follows that there is a set of density Sg at 
least 1 — 5 where {(pUpU^)Y < 5. 

Now we write, from (E2B, 


(7.22) 




{lp'^Uj,Uj)Q - {S(pU^j,u’^j)Y 

{S{1 - (p)u],Vp^Y 


The denominator converges to a positive limit along a density one subsequence S', 
which by (TT^ increases as either ip oi: 1 — (p increases. On the other hand, along 
S'g n Sg, the numerator has oscillation at most 3(5. Hence along the subsequence 
S n S^ n S^', the sequence of norms ||uj||i 2 (y) has oscillation at most a constant, 
depending only on H, times 5. 

By Lemma lT^ this proves that the sequence of norms has a limit along a density 
one subsequence. To compute the limit, we let the support of ip increase to H, and 
the support of (p shrink to E. Then {ip'^Uj,Uj)a 1, {Stpu^ vPj) —> 0, and by 17.31) . 
(TttHI) and rr^ . the denominator of ifT^ converges to 


a;Neu(<5'(l - (p)) ^ 1 

^Neu(Id) CV]SJeu(Id) 


along a density one subsequence. Hence, {u^,u^) converges to a;Neu(Id) along this 
subsequence. This proves (IT^ . and completes the proof of Theorem o for Neu¬ 
mann boundary conditions. 


8. The Robin boundary condition 


The proof for Robin boundary conditions, d,^u = ku, for a real function k on the 
boundary, is almost identical. In this case. Green’s formula gives us 

t N f dGo{x,y;X) i ■ o 

u{x) = J - — - u(y) - Go(x,y,X)—da{y), xG interior H. 

Using the boundary condition we find 

u(x) = J “ K{y)Go{x,y; X)'ju{y) da{y). 


X G interior H. 
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Let Efi be the operator of Proposition KT\ Letting x approach the boundary as 
before, we find that 

u] = F^.u], F^ = Fn + Enon. 

By Proposition^^ the operator Eok has the same structure as F, but is one order 
lower; the FIO part of Eon is order —1, instead of zero. Hence, the principal symbol 
of the FIO part of F” is identical to the case of Neumann boundary conditions. 
Consequently, the result, and the proof, are the same as for Neumann boundary 
conditions. 


9. DiRICHLET BOUNDARY CONDITION 

In this section, we modify the Neumann proof so that it works for Dirichlet 
boundary conditions. 


9.1. Local Weyl law. First, we prove the local Weyl law 
(9.1) ^ ^ WDir(^) 

A,<A 


for boundary values of Dirichlet eigenfunctions. Recall that in this case the bound¬ 
ary values are given by = d,jUj\Y- 

We prove in the same way as in the Neumann case. Reductions (i) - (v) 
work without any modifications. Only the calculation in (vi) in which we consider 
even order differential operators is somewhat different. With Ah — where P 

is a positive differential operator of order 2k as in the Neumann case, we consider 
the slightly different 

OO 

eAit) = 

i=i 


We have 


d , fe+i 


(i) 


eA{t) = 




(F*FuU-) 


trPE{t). 


Hence if we apply PD^Dx' to and restrict to a: = F = 0, we obtain a 

kernel whose trace is eA{t)■ 

Using the expression for the heat kernel, with terms 70 , now given by 

we obtain 


eAit) = 21 e-^^^"+\^\"^\piy,r^)\^edyd^dr^ + 

This is the same expression as (II 2.1 111 (with fc = 0), except that we have an extra 
factor of from the two normal derivatives in t. Following the calculation in the 
appendix, this extra factor has the effect of changing m to m + 1 , and introducing 
the factor I — |(Cp in the integral over |C| < I, changing the Neumann measure into 
the Dirichlet measure. Hence we end up with 



fc+i 


iA{t) t-(t+'=+l)(| + /c)(^ + fc-l) . . 


(-)-- 

2 ( 47 r) 2 vol(S'"' 


B*Y 


a{A) dpBir 
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where ~ means up to an error which is 0{t ("/2+”i+i/2)^ ^s t ^ 0. We conclude 
that 

eA{t) = (47rt)~"/^ UQn-n / t 0. 

V0H;5 ) Jb-Y 

Then, the Karamata Tauberian theorem gives (19.111 . 


9.2. Quantum ergodicity. We now sketch the proof of quantum ergodicity in the 
case of Dirichlet boundary conditions, which is the statement that 

^ WDir(2l), 

along a density one subsequence of integers j. 

A crucial change from the Neumann case is that the formula FhjU^ = in the 
Neumann case gets replaced by 

(9.2) = -u'^^ 

in the Dirichlet case, due to the fact that the normal derivative falls on the oppo¬ 
site coordinate. Hence we now need to use the adjoint of Lemma lb.21 which for 
convenience we state here. 


Lemma 9.1. Suppose Ah is a pseudodifferential operator of order zero, withWF' (A) C 
TZ^. Then, one can express 

{Fhf AhiFf)^ = Bh + Sh, 


where ||>S'ft,||i 2 ^j ^2 < Ch, and Bh is a pseudodifferential operator of order zero with 
symbol 


(9.3) 


a{B) 


7 HlhiP ^{q)), q&Tl ^ 

0, qin-^. 


Hence the Egorov operator T (see section EJ gets replaced by 

Step 1. The first task is to choose a pseudodifferential operator R whose symbol 
is approximately equal to the invariant function c^{q)~^ = c(l — where c 

is as in (trroi) . Given ei > 0, we may choose R so that its symbol is supported in 
B*Y, and such that 


(9.4) 


j\(j{R) - c-i{q) ^|^7(q)dcr(g) < ei. 


We then prove the following analogue of equation GS: 


(9.5) 


lim sup —TTv 
A^oo N{\) 


E YfHAi - 


Xj <c.x 


‘^Dir(A) f, 


Ua 


a;Dir(Id) ^ 


IIl2(v)) 


< Ce, 


where C depends only on H and A. To do this, we define 

-a{l3-\q))^{l3-\q)), 




N ^ 7 ( 9 ) 


(9.6) 
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an averaged version of the symbol in 19.dll . (Here it is not necessary to introduce 
the function a" as in the Neumann proof.) By the mean ergodic theorem, 

^ cuDiriA) in L^. 

In fact, since is uniformly bounded in we also get convergence in for 

2 < p < oo, in particular for p = 6. Since y”^/^ g L^(B*Y), this implies that 

yi/2Q(.;v) ^ in L^- 


Thus, we can choose an integer N such that 

II (a^^^y - ca;Dir(^))y“^^^|| < e. 

Choosing a' as in section [7| we obtain the analogue of (TmUl : 

(9.7) ||y^/^(a'^^^ - cy'^woir)^)) II 2 < 4e. 

The rest of Step 1 then follows the argument in sectional 

Step 2. Again the argument in section 0 adapts with minor changes. Since the 
boundary value here is the normal derivative of the eigenfunction, the value at 
an interior point is given by 

Uj{z) = J^Eh^{z,y)u^j{y)da{y), 

where is the free Green function at energy h~'^, instead of 17.1 fill . Thus, 

, Uj) = {El.'ijj^Ehju] , mJ), 

so we need to analyze the operator E^.i^'^Eh^. Lemma 1711 srets replaced by 

Lemma 9 . 2 . Let (j) be a smooth function on Y which is identically one near E. 
Then 

(i) h~‘^ElfP'Eh{l — (f) is a semiclassical pseudodifferential operator of order zero 
on Y; and 

(ii) For each 1 < i < m, h~^El'ip'^Eh is bounded on L'^lY), uniformly in h. 

(Hi) The limit of ujY)i-[{h~^E*ifAE{1 — </>)), as the support of (f shrinks to E, is 


(9.8) —^ / f^^{z)dz. 

vol(S2) 

The proof follows the same lines as the proof of Lemma 1711 In the proof of (iii), 
the lack of two normal derivatives in the operator is exactly compensated by the 
change in measure from /iNeu = y“^dcr to pDir = ydcr. 

The proof is completed by writing 


(9.9) 


. b ||2 


\L^(Y) 


- Xj '^{S(j)u),u))Y 
\~'^{S{1 - (l))u),uf^Y 


and using Lemma [7?71 and the result of Step 1, as in section [Tj 
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10. IP^-ROBIN BOUNDARY CONDITION 
In this section we consider boundary conditions of the form 

d^u = Ku, 

where K is & self-adjoint pseudodifferential operator of order 1 with nonnegative 
principal symbol on the boundary. We shall denote the operator Hk in this case. 
Recall that it is defined by a quadratic form 
The local Weyl law for Hk is as follows. 

Lemma 10.1. Suppose we impose the boundary condition above. Let k denote the 
principal symbol of K. Let Ah be a semiclassical operator of order zero on Y. Then 


( 10 . 1 ) 


{Ah^u],u)) =ujk[A) = 

\,<\ 


vol(S'"'“i) vol(fl) 


B*Y 


a(g)7(g) 
l{qY + k{qf 


da{q). 


Proof. Identical in structure to the proof of Lemma o 


□ 


Notice that if iL = 0 then this reproduces the result of the previous section. 
Consider now (for a manifold with smooth boundary) K = a Ay . As a —> oo, 
this boundary condition approaches Dirichlet, in some sense. To compare to our 
Dirichlet result, we need to account for the fact that here we are considering the L^ 
norm of the function at the boundary, instead of the normal derivative. To change 
to the normal derivative is equivalent, by our boundary condition, to changing to 
aAy^u, so to compare, we need to multiply the weight factor in (IK).Ill by 
If we do this, then as a ^ oo we get 


lim [ 

a^oo Jk,y 


a{q)i{qW\s? 

'y{q)'^ + a2|,y|2 


dcr{q) = [ 

J B*Y 


a{qh{q)da{q), 


which is the Dirichlet measure. Thus, for this family of boundary conditions the 
corresponding measures LUaK interpolate between the Neumann and Dirichlet mea¬ 
sures. 

In order to adapt the arguments above to the T^-Robin boundary conditions, 
we need an operator, which is essentially a FIO, which leaves the invariant. 


Lemma 10.2. Let e > 0 and N be given. Then there is a subset V of B*Y x 
B*Y such that both projections Vi and V 2 onto the left and right factors of B*Y 
have measure less than e, and a operator , acting on Lf{Y), depending on the 
parameter h > 0, such that 

(i) leaves the invariant: 

= u?, 

(a) there is a decomposition 

F^ = F^ + F^ + F^ 

such that {Ff^)h is a semiclassical FIO of order 0, associated to the billiard map, 
{Ff^)h is a semiclassical pseudodifferential operator of order —1, and {F^)h has 
wavefront set contained in V; 
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(Hi) the symbol of in Cbiiuard \ V is given by 


hiPiq)) 7(g)-»fc(g) 

V 7(g) 7(/3(g))+ *fc(/3(g))’ 

(iv) If Ah is a semiclassical pseudodifferential operator with microsupport disjoint 
from V 2 , then for any —N < m < N, the operator = A^"^^ + S, 

where A^"^^ is a semielassical pseudodifferential operator with symbol 


7(/3™(g)) 


7(g)^ + k{qf 


7(/3’"(g))2 + k{(3'^{q)Y 
and \\Sh\\L^{Y)^L‘^{Y) = 0{h). 


7 ( 9 ) 


■a(/3™(g)) 


Proof. We first prove these assertions in the case = 1. We start from Green’s 
formula: 


u{x) 


dGo{x,y;X) 

-- u{y) - Go[x, y] A)— da{y), 


X € interior Q. 


Using our boundary condition, we find 

u{x) = J ^\ {y) - Go{x,y,X){Ku) da{y). 


Now we let x tend to the boundary. From the normal derivative of Gq we get \5x{y) 
plus half the kernel h = X~^. Using Proposition ^21 given any neighbourhood 
U of the set S x {TZ^)' U (7?.“^)' x S, we can decompose F = Fi + F 2 + F 3 , where 
Fi is an FIO, F 2 is a pseudodifferential operator of order —1 and F 3 has operator 
wavefront set contained in U. 

The other term is more complicated, since we are composing the ‘single layer 
potential’ F, given by the kernel of 2Go restricted to the boundary in both variables, 
with K. It is necessary to introduce a more elaborate decomposition of E in this 
case, since K is not a semiclassical operator, and there is no nice calculus containing 
both homogeneous and semiclassical operators. However, if T is a semiclassical FIO 
whose operator wavefront set is disjoint from the zero section in the right variable, 
that is, if 

WF'(r)c{(g,g')|gVo}, 

then TK is a. semiclassical FIO with symbol t{q,q')k{q'). Here, k is the principal 
symbol of K. This follows from Taylor m, chapter 8, section 7, for example. 

Thus, we let V be any neighbourhood of F U (Z U P{Z)) x Z, where Z is the zero 
section of B*Y. (Notice that both U and V, as well as both their projections, may 
have arbitrarily small measure.) We decompose F = Fi + F 2 + F 3 , where Fi is an 
FIO of order —I microsupported away from the zero section in the right variable, 
F 2 is a pseudodifferential operator of order — I microsupported away from the zero 
section in the right variable, and WF'(F 3 ) C V. 

Thus, we have 


(10.2) = J (Fi + F 2 + F3)u 5 + (Fi + F 2 + E3)Ku^j. 

By Lemmathe symbol of F 2 inside B*Y, and outside V 2 = proj 2 (U), is —* 7 “^. 
Thus, E 2 K is a pseudodifferential operator of order 0 with principal symbol —i^~^k 
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outside V 2 . We rewrite this equation as 

(10.3) (Id -E 2 K^ u) = J {F1+F2 + F 3 )u) + {Fi + F:i)Ku). 

Since the symbol of F2K is imaginary, there is a pseudodifferential operator T such 
that T{ld—F2K) = Id—£^ 4 , where i ?4 is order —1. Moreover, the symbol of T on 
B*Y \ V 2 is (I + Thus, we get 

(10.4) = r(Fi +F2 + F3)u] + T{F^ + F:i)Ku] + F^u], 

Now we let 

F^{x,y)=TFi+TFiK 

(10.5) F^{x, y) = TF2 + Fi and 

F^{x,y) = TF^ + TE^. 

Then we have 

( 10 . 6 ) u] = {F^ + F^ + F^)n,u], h, = Xf, 

where F^ is a semiclassical FIO with principal symbol 


_ I _ I 7 (g) ^fc(g) ^ ^ j liPiq)) -i{q)-ik{q) 

l + ik{/3{q))/j{(3{q))\\^-f{/3{q)) y/j{q)j{f3{q))) \/ l{q) l{(}{q)) + ik{f3{q)) 

on Cbiiiiard \ V; F^ is a semiclassical pseudodifferential operator of order —1, sup¬ 
ported outside V 2 ; and F^ is an error term with microsupport contained in V. 
This proves (i), (ii) and (iii). Consequently, the symbol of {F^)*AF^^, for A mi- 
crosupported within V 2 is 


7(/3(g)) 7(g)^ + fc(g)^ 

lif3{q)y + k{f3{q)y 7 (g) 


a{(3{q)). 


Statement (iv), for = 1, then follows from the calculus of wavefront sets and the 
symbol calculus for FIOs, as in Proposition lb. II 

To prove statement (iv) for arbitrary N we need to enlarge the set P to a set 


y(JV) _ I g ^11 — P^iq") for some 0 < m < TV— 1}. 


The projections onto the first and second factors of this set have small measure if V 
has projections of small measure, since /3 is measure-preserving. Then (iv) follows 
as in the proof of Proposition lti.il □ 


This Lemma allows us to run Step 1 of the argument in section[7|for the ih^-Robin 
boundary condition. 

To run Step 2, we again run into complications caused by composition of homo¬ 
geneous and semiclassical pseudodifferential operators. In the present case, we get 
instead of iTTTTIi 

(10.7) = {{Gh, + Eh^KYij'^{Gh, + Eh^K)u],u'’YY ■ 

Here we need to analyze the operator {Gh + EfiK)*ip^{Gh + EhK). 

Lemma 10.3. Let cj) be a smooth function on Y which is identically one near E, 
and let Q be a semiclassical pseudodifferential operator with microsupport disjoint 
from the zero section. Then 
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(i) {G+EQK)*'ip^(G+EQK){l — 4>) is a semiclassical pseudodifferential operator 
of order zero on Y; and 

(a) For each 1 < i < m, (G + EQK)*'ip'^{G + EQK) is bounded on L^(Y), 
uniformly in h. 

(Hi) The limit o/a;Neu((G + EQK)*ip^(G + EQK){1 — (f)), as the support of (f 
shrinks to E, and the symbol of Q increases to 1, is 


( 10 . 8 ) 


1 

vol(ll) 



dz. 


The proof of this Lemma follows the same lines as that of Lemma o and 
Lemma The novelty here is in dealing with the cutoff operator Q, which 

gives us additional error terms 


{u],G*f:'^E{lA-Q)Ku'>^, [u], {E{ld-Q)Kyif‘^{EQK)u^^, 
and [u], {E{ld-Q)KYy'^{Elld-Q)K)u^y 


and their adjoints. We need to show, with the first term for example, that 

Y. {u).G*y^E{ld-Q)Ku'y = 0 . 

Aj<A 


Writing this as 


lim sup 

A —yoo 


1 


Y {{ld-Q)E*y^Gu],Ku^y) 

\j<\ 


< lim sup 

A —yoo 


1 


Y €-^\\h-\lA-Q)E*y^Gu]f + e\\hKu]f 

Xj <A 


we see that this can be made arbitrarily small by choosing Q suitably, since the hrst 
term is controlled by u!K{(idL—Q){h~^E*)G) which can be made arbitrarily small, 
and the second can be made small by choosing e small. To deal with the third 
term, we need to commute K and {E(ld—Q))*tf^E(ld—Q). This again causes 
difficulties since one is a homogeneous, and one a semiclassical, operator. However, 
Lemma ll2.til allows us to commute the operators up to an error with an o(l) operator 
norm as h ^ 0. This allows us to complete Step 2 for the 'I'^-Robin boundary 
condition, which completes the proof of quantum ergodicity in this case. 


11. NONCONVEX DOMAINS 

Here we adapt the argument to nonconvex domains. The problem with the ar¬ 
gument above for nonconvex domains is that the FlOs Ei^h and Fi^h from Section^] 
have a canonical relation which is larger than the billiard relation, since it relates 
points on the boundary which are connected by a straight line even if the line passes 
outside the domain. In this section, we show that one can manufacture an invariant 
operator F/, which satisfies the conditions of Proposition lO The argument then 
proceeds as before. 

We do this by modifying the metric outside the domain, while keeping it Eu¬ 
clidean inside. Let 6 be a smooth, compactly supported nonnegative function on 
R" which vanishes on H. Consider the metric 


9s — (f T 5^)llEuclidean OR 1^ ■ 
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For sufficiently small s, no geodesics of gg starting at a point in have conjugate 
points in n. Let Gs{z, z, h) = (A^ — z') denote the kernel of the 

outgoing resolvent of the Laplacian on R" with respect to the metric gs ■ 

Proposition 11 . 1 . The kernel Gt{z, z', h) may be represented 

Gt{z, z', h) = z', h), z ^ z'. 

Proof. We obtain the resolvent from the heat kernel by integrating against the wave 
kernel 

/ OO 

Recall that the wave kernel is a family of FIOs, with canonical relation given 

by 


{{z, z', C') I 3 a gs-geodesic 7 of length t with 

z = 7(0), C = C 7 ( 0 ), z' = 7 (s), C' 
It can therefore be written in the form 


-C7(s)}. 


where a is smooth in all variables and symbolic of order —1 in and, Z'{z,z',s) 
is the coordinate of z' in ^g-normal coordinates centred around z. Putting C, = A^, 
we find that 


(A,-(A + 70f)-i =A" 


dt / d^ e 




For |t| > To, smoothness of the wave kernel on fl x fl implies, by repeated integra¬ 
tions by parts in t, that the corresponding kernel is residual. For 0 < to A |f| < 2b, 
where Tq is larger than the ^s-diameter of fl for each s, the integral is a semiclassical 
FIO with canonical relation 


{{z, C, z', C) I 3 5 s-geodesic 7 with z = 7 ( 0 ), C = 7 ( 0 ), z' = 7 ( 5 ), C' = - 7 (s)}- 

This is parametrized by the global phase function dists(z, z'). Finally, for |t| < to, 
where to is sufficiently small, the wave kernel is identical with the free wave 

kernel due to finite propagation speed, so G's(z, z', h) is microlocally identical 

with the free resolvent kernel near the diagonal. This proves the lemma. □ 


To deal with nonconvex domains, we define 

Kiy,y') = ‘^-^Gs{y,y,h), y,y'eY, 

and average over s to obtain the operator 

(11.1) Fn= [\is)Ffds. 

Jo 

Here x is smooth, supported in (0, b), and nonnegative with integral 1. The aver¬ 
aged operator Eh is defined similarly. 

Let us say that a point {]j,g,y',g'), with g = dy\y — y'\, g' = —dy'\y — y'\ is a 
spurious point of WF'(Tb) if the line yy' leaves H. We also recall the definition of 
S in KTh . Then the following analogue of Proposition 14.21 holds. 
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Proposition 11 . 2 . (i) The operator is an invariant operator for the Neumann 
eigenfunctions: 

Fnu]=u). 

Moreover, the analogous statements for other boundary conditions are valid for the 
corresponding averaged invariant operators. 

(a) Let U be any neighbourhood of E x {TZ^)' U x S U x (7?.^)'. 

Then for a suitable averaged invariant operator F^, there is a decomposition 

Fh = Fi,h + F2^h + F^,h, 

where F\ is a Fourier Integral operator of order zero associated with the canonical 
relation Cbuiiard; ^2 is a pseudodifferential operator of order —1 and F^ has operator 
wavefront set contained in U. 

(Hi) The symbol of Fi ^ is the same as for Ff,. 

Proof, (i) Each is an invariant operator since the metrics gt all coincide with 
the Euclidean metric inside LI. Hence any average of the Ff^ is also an invariant 
operator. 

(ii) If 5 > 0 somewhere on yy', then d/ds{dists{y,y'y) > 0 for small s by the 
positivity of b. Hence, in the integral cnj, the phase is locally non-stationary, 
implying that the kernel is locally residual. 

A compactness argument shows that if the support of b is sufficiently close to 
LI, then any spurious vector q G WF'(F;j) is either in U or the line yy' joining the 
endpoints meets the support of b, in which case q ^ WF'(F/i) by the remarks of the 
previous paragraph. 

As for the non-spurious wavefront set, the decomposition follows from the struc¬ 
ture of Gs{z, z', h) given in Proposition lll.il and the proof of Proposition ^21 

(iii) Let {y,ri,y',r]') G WF'(E’?i), such that the line joining y,y' remains inside 

LI. Then the metrics gs are identical in a neighbourhood of this line, hence the 
wave kernels , and therefore, the resolvent kernels Gg are all identical mod¬ 

ulo residual terms near {y,y'). Hence all the Ff are microlocally identical near 
{y,r],y 'which implies (iii). 

□ 

Lemma lb . 21 now holds for a suitable averaged invariant operator Fh. Hence, the 
arguments of Section [3 go through, proving the result in the nonconvex situation. 

12. Appendix. The heat kernel 

We begin by exhibiting the asymptotic expansion of the heat kernel for a domain 
with smooth boundary, first under ^'^-Robin boundary conditions and then under 
Dirichlet boundary conditions. 

12.1. Neumann or ik^-Robin boundary condition. The asymptotic expansion 
of the heat kernel as t 0 may be obtained from the asymptotic expansion of the 
resolvent (A^ — as A —> oo within a sector disjoint from the positive real axis. 
The construction of Seeley in |Sel| , |Se2| goes through for a boundary condition of 
the form 


(12.1) 


d,,u = Ku, KG'I’^{Y), k = ai{K)>Q. 


QUANTUM ERGODICITY OF BOUNDARY VALUES OF EIGENFUNCTIONS 


39 


In local coordinates, the Seeley parametrix takes the form 


( 12 . 2 ) 


i=o 


I e-^^’^e^^y-y'y^d.2-, (X, y, 77, A) d^dy^. 


We obtain the heat kernel e by performing a contour integral: 

e-‘^^ = i X e-‘^(AB - A)-idA, 

where F ‘encloses’ the spectrum of As, eg 

r = {—L + se“*®, s G [ 0 ,00)} U {—L + se*®, s G [ 0 , 00)}, 0 < 6 < tt. 


Substituting the series for i?jv(A) in for the resolvent, we obtain a parametrix for 
the heat kernel of the form in a single coordinate patch 


N 


E ( 2 -)- 


( 12 . 3 ) 


i=o 


-n / i{y-y')-ri 


lj{x,y,^,ri,t)d^dT] 


( 27 r)-" J e^<^y-y'>ye-^^'^Sj (x, y, 77, t) d^dr]\. 


After patching and summing over coordinate patches, we denote this operator 
This is a good approximation to the heat kernel as t —> 0 , in the sense 
that the kernel of the difference, 

( 12 . 4 ) gAr(t, x, y) = {x, y) - ZN{t){x, y), 


satisfies estimates 


(12.5) |D“D^9Ar(t)| < C't-(’^+“+/3)+Af+l)/2e-'5|‘^-!^l"/* 

for some i 5 > 0 . We will be interested mostly in the top terms in the expansion of 
the heat kernel, corresponding to 70 and Sq. To write them in the simplest possible 
way, we assume that we have chosen coordinates so that \dx\ = 1 and dx and dy^ 
are orthogonal at the boundary, and so that the lines {77 = constant} are geodesics 
close to the boundary. Then, the metric gb takes the form 

o\ 

0 

0 ... 1 ] 

where h is the induced metric on the boundary. Let us denote j<n-i by 

\riW- Then the top two coefficients C2 and c?2 in the parametrix for the resolvent 
with boundary conditions (tmil are 


/ 

V- 


( 12 . 6 ) 


C2{x,y,^,r],X) 

d2{x,y,^,r],X) 


ie+Mi-x)-^ 

— k{y, 77) e ^ 
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Recall that k{y, rj) is the principal symbol of RT, hence is nonnegative by assumption, 
and homogeneous of degree one. Hence has no more singularities than it does 
for the Neumann boundary condition k = Q. Essentially for this reason, the Seeley 
parametrix works just as well for 'I'^-Robin boundary conditions as it does for 
Neumann boundary condition. 

The corresponding top terms 70 and 5q for the heat parametrix at a; = 0 are 


(12.7) 




= ^-t{e+\v\l) 

i^- k{y,r]) 


Next, for piecewise smooth domains, we show 


Lemma 12.1. Suppose that C R" is a piecewise smooth domain. Then for 
Robin boundary conditions, the heat kernel admits Gaussian bounds 

for t < 1, 

(12.8) |e-*^^(z,z')| < 


for some constants C,S depending on fl. Moreover, in the regular part Y° of the 
boundary, the local expansion is valid. 


Proof. The operator Ak is the self-adjoint operator associated with the quadratic 
form (lOl . When K is nonnegative, this quadratic form is larger than the cor¬ 
responding form for the Neumann boundary condition {K = 0), so (112.811 follows 
from [D], Theorems 2.4.4 and 3.2.9. 

To prove the second part of the Lemma, we fix a point yg G Y° and extend Q 
to a manifold D H with smooth boundary so that the boundaries of Q and fl' 
coincide on the support of the kernel of K or n, and on a neighbourhood of yo. 
We can then choose a boundary condition B for fl' which agrees with the original 
boundary condition wherever the two boundaries coincide. Since Hq/ has local 
asymptotics, it is enough to show that Hqi — Hq has trivial asymptotics near 1 / 0 ■ 
The heat kernel Ha{t, z, z') for for the boundary condition B can be expressed 
in terms of that for by 

Hait,z,z')=B:Q>{t,z,z')- f [ Ha{t-s,z,y){BHa>{x,y,z'))da{y)ds. 

Jo JY 

Notice that BHQi{x,y, z') = 0 for ?/ near j/q- Hence the Gaussian bounds for Hq 
and Hqi (together with all derivatives of Hqi) imply that the difference between 
the heat kernels has trivial asymptotics, as required. □ 


Our next task is to show that for 4>^-Robin boundary conditions, the operator 
E{t), obtained by restricting the heat kernel to the boundary in both variables, is 
trace class. 

To do this, we prove the bound 
(12.9) < CAJ. 

Note that Ax has no null space since K is positive, and Ax has discrete spectrum 
since the quadratic form is larger than the Neumann quadratic form. Hence Ajf 
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exists, and it is certainly a bounded map from L^(r2) into the form domain 
We conclude that the eigenfunction uj, which satisfies 

satisfies ||wj||//i(n) < Hence by the Sobolev trace theorem, < CXj. 

Since the kernel of E{t) is given by 

( 12 . 10 ) 

3 

and converges for each t > 0, it follows that E{t) is trace class for all 

t > 0. 

Moreover, if H is a differential operator on Y with coefficients supported inside 
y°, of order k, it is easy to show that 

\\Au%.^y^ < CAf 

Since A^“''^ converges for any fc, we have by Lidskii’s theorem that AE{t) 

is trace class, and the trace is given by the integral of the kernel of AE{t) on the 
diagonal. 


Next, we calculate the trace of PE{t), where P is a positive, order 2k differential 
operator on Y with symbol p{y,ri). Decomposing = El^it) + q^it) as in 

(tmi . we first substitute for E{t) the operator obtained by restricting the kernel 
of Hq to the boundary in both variables. This we can compute directly: we get 

(12.11) (271)"”' f - -dyd^dy. 

To compute this we make a linear change of variable from rj to Q = Q{y, rf) which 
is an orthonormal basis for h^^{y) for every y. Write /x = det(/i*l)“^/^, so that the 
Riemannian measure on T is ^ dy. Under this transformation, p{y, rf) changes to 
P2k{y,0 +p2/c-i(y,C)i say, where p 2 k is homogeneous of degree 2k in ( and P 2 k-i 
is a polynomial of degree 2fc — 1 in (f. We obtain 

(27r)"" J {p 2 kiy, C) + hk-iiy, o) dy d^ dC 

= ( 27 r)"”y +P2k-iiy,0) ^^2 ^^^2 k-dydS,d(f. 

The —2if,k term can be dropped because this gives an odd integral in If we then 
change to polar coordinates, then we get 

(12.12) {2n)- J" dr J ydy J p,,{yX)^^ + ... 

° l(C.C)l=i 


Here the dots represent the contribution of P 2 k-i- H we integrate by parts k times, 
we find that 


— tr^n—l-\-2k 


dr = {n/2 + k — l)(n/2 + fc — 2 )... {n/2)t 


-k 


e-*" r^-^dr. 


Turning this back into an integral over R" gives 

1 


g-tr- ^n-l 


dr = 


VOl(S'" 1) jRn 


/ 

J'MJ' 


dx = 


vol(S'"-i) 


^nj2^-n/2^ 
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Writing = 1 — |Cp on the unit sphere and representing the measure as (1 — 
we get 


(f + fc - l)(f + fc - 2)... 

(47r)’i/2 



ICI<1 


P 2 k{yX) 


4v^r^ 

l-|CP+fc 2 


dc + • • • 


The extra factor of 2 comes from the fact that there are two values i\/l — |CP on 
the unit sphere for each ( with |(^| < 1. It is now clear that the contribution of the 
omitted term is Thus, we have shown that 


(12.13) 
ti PE{t) 


n n At 2 

^2^ ~ ^■■■^2^(47r)"/2vol(5'"-i) 


aiP)dtiB +0(t-t-'=+^) 


B*Y 


as t —> 0. 

This result holds in particular for multiplication operators cj), supported away 
from S: 

(12.14) tr(^g(t)^ (4^)n/2 / (l>da + 0{t-^ + ^). 

Y 

The Gaussian bound iinii shows that this holds for all multiplication operators 
(/), not necessarily supported away from the singular set. Hence (112.1311 holds also 
for P equal to the identity operator: 

2f-"/2vnlWi 

(12.15) tr Fit) = ^ t 0. 


12.2. Dirichlet boundary condition. The Seeley parametrix for the Dirichlet 
boundary condition on a manifold with smooth boundary still takes the form (tTOl . 
but the top terms 70 and 60 now are 

(12.16) ,, 

Mx,y,^,V,t) = +l'/lh)e-*€^. 

Lemma imi gets replaced by 

Lemma 12.2. Suppose that H C K" is a piecewise smooth manifold. Then the local 
expansion of the Dirichlet heat kernel {z, z') is valid on any compact subset 

of the regular part Y° of the boundary. 

Proof. This is proved as for Lemma nm using the bound 

(12.17) |e-*^“'(z,z')| < (4^t)-"/2e-l"-"'lV4t_ 

which follows by comparison with the heat kernel on R". □ 


Next, we prove an identity for Dirichlet eigenfunctions on Lipschitz domains 
which is elementary for manifolds with smooth boundary. 

Lemma 12.3. The Dirichlet eigenfunctions Uj satisfy the equation 

(12.18) ~ J [x ■ r')\Uj\^ da. 
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Proof. For a smooth domain this formula was first proved by Rellich; see | HT| for 
a proof. For a Lipschitz domain, the conclusion is still valid; see |Ve| . sections 1 
and 2. □ 


A similar proof, with the vector field X replaced by one which is pointing out of 
the domain everywhere, shows that 

||u5f < CXl 

It follows from this that the kernel E{t) is trace class, the argument being the same 
as in the Neumann case. 

Next we show 


Lemma 12 . 4 . Let P be either the identity operator on Y , or a differential operator 
of order k on Y supported away from E, with prineipal symbol p{y,rf). Then the 
trace of PE(t) satisfies 
(12.19) 

trPEm = (f + + k-l)... + o(r ?-*-■). 

B*Y 

Proof. We first consider the case where the coefficients of P are supported away 
from S. We compute as for the Neumann case, using the expression (I12.dll for 
the heat kernel, with terms 70 , So now given by cna, and remembering to take 
normal derivatives in both variables. We obtain 

tr PE(t) = {2 tt)~’^ J +P 2 k-i{y,C))‘^^'^ h-dyd^dr]. 

Following the Neumann calculation, the P 2 k term gives irroTi . while the P 2 k-i 
term gives an error. Note that the effect of the Dirichlet boundary 

condition is to change fc to k + 1, and to introduce the factor into the integral 
above. This factor becomes 1 — in the integral over |C| < 1, changing the 
Neumann measure into the Dirichlet measure. 

The result for the identity operator is a special case of the following lemma. □ 


Lemma 12.5. For every measurable set A C Y, 

1 




1V(A)'P 


j, 2 . 2vol(A) 

u f da = - 

nvol{il) 


Remark. This result was proved for domains with smooth boundary in 
Proof. By the Karamata Tauberian theorem, it is sufficient to show that 


E< 




.b\2 . 


2vol(A) ^_„/2 


IA 


da t t 

(47r)"/2n 


0. 


By integrating in t down from t = 1, it is sufficient to prove that 
..-tA? f L,b|2T_ .-n/ 2-1 


E' 


da • 


(Itt)”/' 




0 . 


This is the same thing as asking whether 
(12.20) / E{t,y,y)da{y) 


vol(A) ^_„/ 2 _i 


(47r)"/2 


t^O. 
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Let US consider the quantity 

(12.21) f {x-iy)E{t,y,y)da{y). 

J A 

By (I12.19II . if A is compactly contained in Y°, then 112.2111 has asymptotics 
(4.)ivoi(5-4 J + 

B*A 

Since ^ ^ 

[ v^l — dr = — [ r"~^ dr, 

./o ^0 Vl — 

WDir(l) = a;Neu(l)/?T-, SO 112.2211 is given by 


(12.23) (47r) 2 ^ ^ {x ■ v) da + o{t = ^ i). 

Ja 

Next consider the case when A = Y. Then using 112.1811 . 112.2111 is equal to 

(12.24) J {x • v)\u^jf da = 2^^ X^{uj,Uj) = 2 ti Ae~*^. 


This has asymptotics 

(12.25) 2trAe-*^ = -24tre-*^ = t-/2-i + ^(t- 

dt (47r)"/^ 

= {4Tr)-^(^J X-+o(t-t-'=-i), t ^ 0. 

Now suppose for a moment that Q is starshaped. Then we can choose an origin 
so that X ■ v is positive on Y, and this together with the positivity of E{t,y,y), 
the fact that Y° has full measure, and implies that ilTT^ holds for all 

measurable A. 

Now we note that property 112.2011 depends only the asymptotics of E{t,y,y) 
as t ^ 0. This is local in y, so (ITT^ therefore holds for locally starshaped 
domains. But this includes all Lipschitz domains, so we have finished the proof of 
the lemma. □ 

12.3. Commutators of pseudodifferential operators. In this subsection we 
prove a Lemma needed in section cni 

Lemma 12 . 6 . Let K he a homogeneous pseudodifferential operator of order m and 
Qh a semiclassical pseudodifferential operator on Y of order 0, both supported away 
from S in both variables. Assume also that Qh has compact microsupport. Then 
the operator norm of [K, Qh] is o{h~”^) as h ^ 0. 

Proof. We prove this in the case of operators on a manifold without boundary. The 
result above follows by localization to a neighbourhood of the support of K. Thus, 
we assume now that E is empty. Also, we restrict to the case to = 0, since the 
general case follows with minor modifications to the argument. 

Let A be some positive elliptic differential operator of order 2 on Y. We first 
prove the result when Qh = 4 >{h?^) for some smooth function (j) with compact 
support. One can expect the theorem to be easier in this case since Qh is a semi¬ 
classical operator, but closely related to the homogeneous operator A which has 
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well-behaved commutators with K. We use the commutator formula and an almost 
analytic extension $ of </) to write the commutator 

[K,cj^{h^A)] 

1 


= [K,h^A](l)'{h^A) + 


2tt 


d^z)ih^A - z)-^[[K, h^A],h^A]{h^A - z)-^ dzdz. 


which follows from the Helffer-Sjostrand formula 

(l){h^A) = — [ d<^>{z){h^A- z)-^dzdz. 

Stt Jc 

The first term may be written 

h^[K, A](l h^A)-^/^{l + h^Af/'^(t)'{h^A) 

(12.26) = h^[K, A]A-i/2ai/2(i + h'^ 

= h{[K, A]A-i/2) {hA^/‘^(l + A)-^/‘^'^^{h^A). 

The first factor in large parentheses is a pseudodifferential operator of order zero, 
the second is an operator with operator norm bounded by 1 and the third is a 
semiclassical operator of order zero. Hence the operator norm of the (112.2611 is 
0{h). 

To deal with the integral over C, we insert the factor AA ^ = {A — h '^z + 
h~‘^z)A~^, and get 

- / M{z){h^A-z)-^AA-\[K,h^Alh^A]{h^A-z)-'^dzdz 

Jc 

= Ch^ j d^{z) ^ Id +z{hJ'A — (hJ A — z)~'^ dzdz 

Recall that |9$(z)| < Cn\ Im2;|“^“^(z)“^ for any N. Hence we can estimate the 
operator norm of this integral by 

Ch^ f d<^{z){l + {z)\lmz\-^)\lmz\-^ dzdz < Ch^{z)-^ < Ch^. 

Jc 

This proves the Lemma when Qh = 4>{h?A). 

To prove the Lemma in general, we choose a function (j) so that <(>(|?7p) is identi¬ 
cally 1 on the microsupport of Qh- This is possible since Qh has compact microsup¬ 
port by hypothesis. Then, for any e > 0, we can find an /iq > 0 and a differential 
operator P, or order m say, such that the operator norm 

WQh - h™P(^(h^A)||i2(Y)^L2(y) < e for all h < ho. 

This follows from the density of polynomials in the C°° topology on compact sub¬ 
sets, as in the proof of Lemma o in Section O Hence, it is enough to prove the 
result for all operators of the form h"‘P(/)(/i^A). (This is rather similar in spirit to 
the usual proof of the Riemann-Lebesgue lemma.) 

To do this, we write 

[K,h'^P(l}{h'^A)] = h'^P[K,(l}{h^A)] + [K,h'^P](l){h^A), 
and consider each term. The second term is of the form 

(12.27) h^A(j){hJA), where H is a pseudo of order m — 1, 
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which can be treated as in (I12.2HII . The hrst we expand as 
h^P[K, 

= (P(l + A)-™/ 2 ^ ((/i^ + /i 2 A)’”/ 2 (l + (^(1 + /i 2 A)^[A:, ./^(/i^A)]). 

The first two factors in large parentheses are bounded operators uniformly in h, 
provided N > ml2. For the remaining factor we write 

(1 + h'^A)^[K, 0(/i"A)] = [K, (1 + /i2A)^0(/i2A)] - [K, (1 + A)^]<^{h'^A) 

and note the the first term has operator norm 0{h) by the first part of the proof, 
while the second is a sum of terms of the form (1TT77II . This completes the proof 
of the Lemma. □ 
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